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Abstract 

The Exact Parity Model (EPM) is, in part, a theory of neutrino mass and 
mixing that can solve the atmospheric, solar and LSND anomalies. The cen- 
tral feature of the neutrino sector is three pairs of maximally mixed ordinary 
and mirror neutrinos. It has been shown that ordinary-mirror neutrino os- 
cillations can generate large neutrino asymmetries in the epoch of the early 
universe immediately prior to Big Bang Nucleosynthesis (BBN). The large 
neutrino asymmetries generically suppress the production of mirror neutri- 
nos, and a sufficiently large Ue asymmetry can directly affect light element 
synthesis through nuclear reaction rates. In this paper we present a detailed 
calculation of neutrino asymmetry evolution driven by the six-flavour EPM 
neutrino sector, focusing on implications for BBN. 
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I. INTRODUCTION 



It has been known for a long time, but not widely appreciated, that parity can be a 
symmetry of nature if the particle content is doubled In this circumstance, for each 

ordinary particle there is a mirror particle of exactly the same mass as the corresponding 
ordinary particle. The set of mirror particles interact with each other in exactly the same 
way that ordinary particles interact with themselves. The mirror particles are not copi- 
ously produced in any laboratory experiments because they either do not couple, or couple 
extremely weakly, to the ordinary particles. In the modern language of gauge theory, the 
mirror particles are all gauge singlets under the standard model G = SU(3)c®SU(2)i(S)U(l)y 
gauge interactions. Instead, the mirror particles interact with a set of mirror gauge particles. 
This is mathematically described by a doubled gauge symmetry of the theory, that is G is 
extended to G (S> G. (The ordinary particles are of course singlets under the mirror gauge 
symmetry.) Parity is conserved because the mirror particles experience V + A mirror weak 
interactions instead of the usual V — A weak interactions. 

The ordinary and mirror sectors can interact with each other in a number of ways. All of 
these interactions apart from gravity can be controlled by an a priori arbitrary parameter. 
Apart from the irremovable gravitational interaction, there are three other ways in which 
ordinary and mirror particles can interact with each other. Two of these are: photon - mirror 
photon (and Z - mirror Z) kinetic mixing 0,^] and Higgs - mirror Higgs mass mixing 0. 
If one demands that the reasonably successful Big Bang Nucleosynthesis (BBN) predictions 
not be greatly disturbed then it seems unlikely that these interactions can have observable 
laboratory implications. See Refs. [^,0 for details.^ 

Neutrinos provide a third possible interaction between the ordinary and mirror sectors. 
// neutrinos have mass then mass mixing between ordinary and mirror neutrinos is possi- 
ble. This leads to very important experimental tests of the exact parity idea. We call the 
G(g)G extension of the Standard Model the Exact Parity Model (EPM) ||,|]. It is, in part, 
an explicit theory of neutrino mass and mixing. It is a candidate for the Standard Model 
extension called for by the solar, atmospheric and accelerator neutrino experiments that 
strongly suggest the existence of neutrino oscillations. Ongoing and future terrestrial ex- 
periments, such as SuperKamiokande, the Sudbury Neutrino Observatory (SNO), Borexino, 
the long and short baseline neutrino oscillation searches and other experiments, will in the 
next few years provide important new clues in the search for a theory of neutrino mass and 
mixing, and will further test the proposed EPM resolution of all of the anomalies. 

Of course, if mirror matter exists then there should be dramatic implications for as- 
trophysics and cosmology as well as particle physics. Some studies suggest that mirror 
matter is an interesting candidate for dark matter. In fact, there is some evidence that 
mirror stars may have already been discovered in the MACHO experiments p|. Another 



^In particular extensions of the exact parity idea, neutral gauge boson kinetic mixing and/or 
Higgs mixing may not be controlled by an independent arbitrary parameter. For instance, in exact 
parity extensions of grand unified models such as SU(5) the kinetic mixing parameter is calculable 
as a function of other parameters in the theory Q . 
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exciting possibility is that Gamma Ray Bursts may be due to collapsing or merging mirror 
stars Of direct relevance to the present paper, however, is the observation that early 
universe cosmology, through Big Bang Nucleosynthesis, structure formation and in the near 
future through detailed cosmic microwave background measurements, should provide im- 
portant new information about the cosmological role of neutrino physics. This information 
may thus also provide a test of the EPM. Indeed, the purpose of this paper is to perform 
a detailed study of the early universe cosmology of the EPM, with particular emphasis on 
BBN. Before embarking on this analysis, we will briefly review why the EPM supplies an 
interesting theory of neutrino mass and mixing, and therefore why some effort to study its 
early universe cosmology is justified. 

It was pointed out several years ago [§] that the EPM provides an interesting theory 
of neutrino mass for one simple reason: the exact parity symmetry between the ordinary 
and mirror sectors forces an ordinary neutrino I'a (a = e, /i, r) to be maximally mixedQ 
with its mirror partner u'^. It is certainly very interesting that the atmospheric neutrino 
observations of SuperKamiokande and other experiments |12| point to the muon-neutrino 
being maximally mixed with another flavour Ux. It is known that cannot be the z/g ||13|| , 



which leaves Ux = and Ux = Vs as the viable possibilities [14,15|, where the subscript s 



denotes a sterile neutrino. The EPM provides a natural candidate for z/^,, namely the mirror 
muon-neutrino v'^. As far as terrestrial experiments are concerned, the v'^ is a sterile flavour. 

Since the confirmation of atmospheric v^^ disappearence by SuperKamiokande, a signif- 
icant amount of theoretical effort has gone into trying to explain the large mixing angle 
observed. This work has focussed almost entirely on the Vx = Vt possibility. It is interesting 
to note that in the immediate past, small interfamily mixing, as observed for the quark sec- 
tor, was considered to be natural also for the lepton sector. With the advent, in particular, 
of the beautiful SuperKamiokande results, this theoretical prejudice is now being criticised. 
Our proposal is completely different from any current effort to realise large — Vr mixing. 
We simply argue that the connection between exact parity symmetry and ordinary-mirror 
neutrino maximal mixing is an especially elegant and simple explanation of the large mixing 
angle observed in the atmospheric neutrino experiments. Since this points to Ux being 



rather than u^-, a neutral current atmospheric neutrino measurement is vital |16 



Intriguingly, there is independent experimental evidence for large angle neutrino oscilla- 
tions from another set of measurements: maximal mixing between the electron-neutrino and 
some other flavour is well motivated by the solar neutrino problem In the EPM, 



the "other flavour" is of course the mirror electron-neutrino fg. Such a scenario leads to 



^Of course this result only holds if the parity symmetry is not broken by the vacuum. In Ref. Q it 
was shown that this occurs for a large range of parameters with just the minimal Higgs sector of one 
Higgs doublet and one mirror Higgs doublet. It was explained in Ref. |1C] that if additional Higgs 
scalars exist, then the parity symmetry can be spontaneously broken with the mirror electroweak 
symmetry breaking scale left as a free parameter. It was argued in Ref. that such a scenario 
could be motivated by the neutrino anomalies. Of course the implications for neutrino experiments 
and early universe cosmology of the model in Ref. [10| are quite different from the minimal case 
considered in the present paper, where the parity symmetry is not broken by the vacuum. 
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an energy independent 50% solar v^. day-time flux reduction^ for a squared mass difference 
greater than about 3 x 10~^° eV^, and to a "just-so" picture for a squared mass difference in 
the approximate range,P|few x 10~^^ — 3 x 10~^° eV^. The most recent solar neutrino data, 



compared with the most recent solar model calculations show that four out of the five 



solar neutrino experiments observe close to a 50% fiux deficit. (The Chlorine experiment 
sees a greater than 50% deficit.) The detailed implications of the solar neutrino situation, 
though surely indicative of oscillations, is not at present as clear as the atmospheric 
neutrino situation. For various reasons, more experiments are needed: the relatively low 
Chlorine result needs to be checked by another experiment, and the cause/existence of the 
apparent distortion of the "boron" neutrino energy spectrum requires further investigation. 
Of particular relevance for the EPM, we will presumably soon find out from SNO whether 
the solar neutrino fiux contains a significant sterile component. 

Notice that no mention was made of the Liquid Scintillator Neutrino Detector (LSND) 
observations in advocating the existence of what are essentially light sterile neutrinos. 
It has become commonplace to motivate light sterile neutrinos from the inability of three- 
flavour oscillations to simultaneously resolve the atmospheric, solar and LSND anomalies. 
We have used this argument ourselves. We would, however, like to emphasise that our 
obsession with the EPM arose from the maximal mixing feature, long before the advent of 
LSND. 

Let us turn, then, to the cosmological implications of mirror neutrinos. The distinction 
between mirror neutrinos and strictly sterile neutrinos, which is totally unimportant for 
terrestrial, atmospheric and solar neutrinos, is of some significance in the early universe. 
This issue will be discussed in depth in later sections. For the purposes of these introductory 
remarks, however, the distinction need not be made. Many of the qualitative features of 
sterile neutrino early universe cosmology pertain also to mirror neutrinos. 

In recent years, the physics of active-sterile neutrino oscillations during and before the 
BBN epoch has been re-examined [E5|-|2n]. Prior to this re-analysis, it had been concluded 



that light sterile neutrinos were cosmologically disfavoured for much of parameter space 
||30| . Focusing on the — > solution to the atmospheric neutrino problem by way of 
concrete example, it had been concluded that the oscillation parameters required would 
lead to the being thermally equilibrated prior to BBN, thus increasing the expansion 
rate of the universe and worsening agreement between theory and primordial light element 



In a very interesting recent paper, Guth, Randall and Serna [19| have pointed out that an 
energy-dependent day-night effect in general exists for solar neutrinos even if the vacuum mixing 
is maximal, thus correcting a misconception shared by the present authors and some of the rest of 
the community. It is not correct to conclude that maximal oscillations out of the "just-so" regime 
always leads to a completely energy independent suppression, because the night-time rate is in 
general energy-dependent due to matter effects in the Earth. 

^There is also an interesting 'window' around brr? 5 x 10"^'^ eV^ which leads to an 



approximate energy integrated flux reduction of 50% and can also explain the distortion of high 



energy ^ ~ 13 MeV boron neutrinos suggested by recent SuperKamiokande data pl[| . 
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abundance measurements. However, it was subsequently realised that the explosive 
production of large neutrino-antineutrino asymmetries or chemical potentials by the active- 
sterile oscillations themselves had not been properly taken into account in the early studies. 
Large neutrino asymmetries generically suppress active-sterile oscillations by making the 



effective mixing angle in matter very small |31]. Detailed numerical work has shown that, 



for a large region of parameter space, the generation through the oscillations themselves of 
large neutrino asymmetries suppresses the production of sterile neutrinos sufficiently for the 
expansion rate of the universe during BBN to be essentially unaffected [p5| , p^ . 

Furthermore, unless the mixing between the electron-neutrino and the other neutrinos is 



really tiny, one expects an asymmetry to develop for z/g's [^. This has a direct effect on the 
rates of the weak interactions processes Ugn e~p and u^p — > e+n which help to determine 
the neutron to proton ratio during BBN. A detailed calculation within a particular neutrino 
mass and mixing scenario is required to work out the magnitude of this effect. It has been 
shown for two different "3 active plus 1 sterile neutrino models" that the generation of a z/g 
asymmetry can be important ||26|j27[] . In a Sec.V, we will for the first time explore this effect 
for mirror rather than strictly sterile neutrinos. 



Aspects of the early universe cosmology of mirror neutrinos were discussed in Ref. p2 
The present paper improves and extends this analysis. Reference |]32| focussed solely on what 
we can call "high temperature neutrino asymmetry evolution". (We will explain precisely 
what we mean by this designation later on.) It showed that the z/^ u'^ and z/e z/g 
solutions to the atmospheric and solar neutrino problems, respectively, were compatible 
with BBN for a large range of parameters. In Ref. , the calculations were carried out in 
the "static approximation" . In the present work we improve on these calculations by using 
the full quantum kinetic equations, rather than the above approximation. In addition, we 
also analyse the "low temperature neutrino asymmetry evolution" that occurs immediately 
prior and during BBN. The size and evolution of the z/g asymmetry will be the main issue 
here. 

Finally, let us remark that the neutrino phenomenology of the EPM is very similar 
to some models employing pseudo-Dirac neutrinos [33|. Many of the implications for early 
universe cosmology will be qualitatively similar to the EPM. There are of course quantitative 
differences because the mirror weak interactions play an important role in the early universe 
through their impact on the matter potential and also because they affect the momentum 
distribution of the mirror neutrinos. We focus on the mirror neutrino scenario in this paper 
because it is arguably much more elegant from a model building point of view. (For example, 
the see-saw mechanism can be invoked to understand the smallness of both the neutrino and 
mirror neutrino masses |^.) It is also theoretically very well motivated because it restores 
parity as an unbroken symmetry of nature. 

The outline of this paper is as follows: In Sec. II we define and motivate the neutrino mass 
and mixing parameters which we will use in our subsequent analysis. We also very briefiy 
review the neutrino asymmetry amplification phenomenon. In Sec. Ill the quantum kinetic 
equations for ordinary - mirror neutrino oscillations are defined and discussed. In Sec. IV 
we compute the region of parameter space where the Ur ^ ^'^ oscillations generate L^^ and 
L,ji^ asymmetries in such a way that the maximal z/^ ^ u'^ oscillations cannot significantly 
populate the mirror i/'^ states. The main issue here is whether or not the z/^ u'^ oscillations 
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can produce compensating L^^ and L^i^ asymmetries such that the matter term for z/^ ^ z/^ 
oscillations becomes unimportant. In Sec.V the low temperature evolution of the neutrino 
asymmetries is studied in detail. The main issue here is the effect of the oscillations on 
BBN. In Sees. VI and VII we comment on the implications of the EPM for the hot plus cold 
dark matter scenario and the anisotropy of the cosmic microwave background. Section VIII 
is a conclusion. 



II. OVERVIEW AND ORIENTATION 



The analysis of neutrino oscillations in the early universe is complicated. In order to 
avoid the pedagogical danger of becoming mired in the full technical detail, we present first 
a short overview. 

There are six light neutrino flavours in the Exact Parity Model: the three ordinary 
neutrinos i^e,^i,T and their mirror partners fg/xr) respectively. In the absence of interfamily 
mixing, the most general neutrino mass matrix consistent with parity symmetry for each 
generation is contained in P[ 





( mi 














\ m2 


ml 1 





+ H.C. 



(1) 



We have assumed Majorana masses for definiteness and simplicity, and one should note that 
the parity symmetry interchanges vl with '^^v'^. The quantity m2 must be real, while mi 
may be complex. However, the phase of mi can, without loss of generality, be absorbed by 
the neutrino and mirror neutrino fields. In the phase redefined basis, the mass matrix is 
diagonalised by the orthogonal transformation. 



J R 




R 



(2) 



We see that the mass eigenstates {v±) are maximal combinations of the weak eigenstates 
(and vice- versa). Obviously it follows that if, as in the quark sector, the mixing between the 
generations is nonzero but small, then each pair of weak eigenstates. 



{Ve,y'e), (j^r^K) 



(3) 



is approximately given by an orthogonal pair of maximal mixtures of the appropriate pair 
of mass eigenstates. We use the notation 



(4) 



for the mass eigenstates. The subscript in the above equation is used to indicate the pair of 
states which relate to the corresponding weak eigenstates. In the limit of no mixing between 
the generations. 
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Of course the exact expressions for Va and z/^ (a = e, /x, r) will in general be a linear 
combination of all possible mass eigenstates when mixing between generations exists. This 
means that all possible oscillations modes amongst the six neutrino flavours are in general 
expected to occur. The assumption of small mixing between the generations, together 
with the necessarily maximal mixing between the ordinary and mirror neutrinos of a given 
generation, implies that intergenerational modes such as v-r ^ v'^ or will have much 

smaller amplitudes than the ^ ^'a modes (in vacuum). The analysis to follow will only 
consider the region of parameter space where vacuum mixing between generations is small. 

In order to proceed, we also have to make a guess about the pattern of mass eigenvalues. 
We will suppose that the neutrino sector is qualitatively identical to the quark and charged- 
lepton sectors, with the masses displaying the standard hierarchy. We will further assume, 
most of the time, that the mass splitting between the parity partners within a given family is 
smaller than the interfamily mass splitting. Putting this together, we have the mass pattern 

- > - m^^. > - m^,_ ■ (6) 

The LSND result suggests that the e — fi mass splittings are of the order of an eV or so, 
although we will also consider smaller mass splittings. If the e — /i mass difference is of the 
order of an eV, then to maintain the assumed mass hierarchy the Ur and i>!j. masses should 
be larger than or about a few eV. A mass in the few eV range would of course make Ur a 
hot dark matter particle. Cosmological closure puts an upper bound of about 40 eV on m^^. 
Analogy with the quark sector suggests that neutrinos in adjacent families, e — fi and fi — r, 
should mix more strongly than e — r. Furthermore, one might guess that a — (3 /a' — (3' 
mixing should be stronger than a' — (3/a — (3' mixing if one believes that the more "closely 
related" are the neutrinos the more strongly they should mix. (Also observe that the parity 
symmetry forces the a — {3 and a' — (3' mixing angles to be equal, similarly the a' — (3 and 
a — (3' mixing angles.) Putting these guesses together with the i/^ — v'e solution to the solar 
neutrino problem and the v'^ solution to the atmospheric neutrino problem, we arrive 

at the parameter space region^ 

^u,- < 1 eV, 10-" eV^ ~ \5ml^,\ = \ml^^ - ml^_ \ ~ lO^^ eV^, 
^^^^+5 'f^Uf,- ~ few eV, 10"^ eV^ ~ = ~ "^I^A ~ 10"^ eV^, 

few eV ~ m^^_^, m^^_ ~ 40 eV, \5ml^,\ = \ml^^ - ml^_ \ < 1 eV^ (7) 

with a mixing angle pattern as described above. 

We wish to calculate the effect on early universe cosmology of neutrino oscillations within 
the EPM. A full six-flavour analysis is a daunting task, even with the parameter space restric- 
tions discussed above. Fortunately, the physics of the problem allows some simplifications to 
be made without sacrificing too much in the way of rigor. In particular, we can build on what 
we already know about the early universe cosmology of active-sterile neutrino oscillations. 

It is useful to start by identifying four qualitatively different epochs: 



^If 3 X 10 ^ ~ |mgg, l/eV^ ~ 10 then the electron neutrino oscillations can have potentially 



observable effects for atmospheric neutrinos, see Ref. [34| for details. 
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1. the Quantum Zeno Epoch, where neutrino oscillations are completely damped; 

2. the High- Temperature Epoch, where large neutrino asymmetries are initially gener- 
ated; 

3. the Low- Temperature Epoch, where decoherence can be neglected; and, 

4. the Big Bang Nucleosynthesis Epoch, where neutrino oscillations impact on light ele- 
ment synthesis. 

We now very briefly, and qualitatively, discuss these epochs in turn. The mathematics 
needed to fully explain this cosmological history is available in previous publications and in 
later sections of this paper. 



A. The Quantum Zeno Epoch 



Neutrino oscillations in the early universe are always to some extent damped through 
collisions with the background medium. As we look back toward the Big Bang, the collision 
rate increases as (below the electroweak phase transition). At sufficiently high tempera- 
tures, collisions occur so frequently that the quantally coherent oscillatory behaviour cannot 
develop. The neutrino ensemble is frozen with respect to its flavour content (Quantum Zeno 
Effect). In addition, the finite temperature contributions to the effective matter potentials 
for many of the oscillation modes are high enough to render the associated matter mix- 
ing angles extremely small. So even with collisions artificially switched off, many of the 
oscillation modes would have tiny amplitudes. 



B. The High- Temperature Epoch 

As the temperature decreases, collisional damping is reduced, and partially incoherent 
evolution of the neutrino ensemble begins. For simplicity, we will in this and the next 
subsection very briefly review the evolution of the a-like lepton number in the somewhat 
artiflcial case where only the z/^ ^ mode is operative. It was shown in Ref. that 
under the influence of this modeQ the a-like lepton number L^^ evolves as per 

- C fL„. + (8) 



dt V " 2 

The a-like lepton number is deflned by 



^ ^'^^ (9) 



^For the purposes of this introductory discussion, the distinction between mirror and sterile neu- 
trinos will often be neglected. 
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and is synonymously called the "a-like neutrino asymmetry" . The quantity rii is the number 
density for species i. Equation (|) holds provided that (i) the squared mass difference Sml^g 
between the neutrinos obeys |5m^g| ~ 10"'^ eV^, and (ii) L^,^ is small. The quantity t] is set 
by the relic nucleon number densities and is expected to be small: ?7/2 ~ 10"^''. The term 
C is a function of time t (or equivalently temperature T). At high temperature it turns 
out that C is negative, so that (L^^ + i]/2) ^ is an approximate fixed point. However 
if Srri^g < [our convention is defined in Eq.([l5|) below], then C changes sign at a 
particular temperature T = Tc, estimated to be |^ 



T,^ie(=^<l^y MeV. (10) 

At this temperature, rapid exponential growth of neutrino asymmetry occurs, unless sin^ 29as 
is very tiny [see Eq.(|lT]) below]. The generation of neutrino asymmetry occurs because the 
t'a — i^s oscillation probability is different from the z/^ — Vg oscillation probability due 
to the matter effects in a CP asymmetric background. As the asymmetry is created, the 
background becomes more CP asymmetric because the neutrino asymmetries contribute 
to the CP asymmetry of the background. This leads to a period of runaway exponential 
growth of the neutrino asymmetry for a large range of parameters p5|j3^n summarised by: 

Sml^ < with \6ml^\ ~ 10"^ eV^, 



10-1° ~ sin^ 29^s ~ few x lO^^ 
10-1° ~ sin^ 26,,, ~ few x 10"^ 



eV2 



l<5mLI 
eV2 



for ordinary - sterile oscillations, 
for ordinary - mirror oscillations. (11) 



(The upper bound in the above equation comes from a constraint on the effective number 
of neutrino flavours, N^^eS, during BBN. We have used N^^^s — 3 ~ 0.6 in this equation 
for illustrative purposes.) We want to emphasise and to state very clearly the following 
fact: Provided the oscillation parameters are in the large range given in Eg. (PUj, the or- 
dinary - sterile (or mirror) neutrino oscillations will generate, at the temperature T^, a 
significant neutrino asymmetry (or chemical potential) from the tiny seed CP asymmetry of 
the background plasma. There is no choice about this, a point sometimes misunderstood in 
the literature: the large neutrino asymmetry will inevitably be generated. Once generated, 
neutrino asymmetries in turn contribute to the effective matter potentials and generically 
suppress oscillations by inducing small effective mixing angles. For typical oscillation pa- 
rameter values within our scenario, the explosive neutrino asymmetry growth begins while 
collisions still dominate the evolution (though they now do not completely damp the oscilla- 
tions). Note that the evolution of lepton number for T < is approximately independent of 
the initial neutrino asymmetries provided that they are not too big (that is, less than about 
10^^). This is because of the approximate fixed point structure which sees L^^ —v/'^ fo^ 
T>Te. 



^In the region of parameter space where |(5m^g| <C 10 ^ eV^, the evolution of the neutrino asym- 
metry is dominated by oscillations between collisions and the lepton number tends to be oscillatory 
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C. The Low- Temperature Epoch 



While neutrino asymmetries develop and evolve, the collision rate continues to decrease in 
a fashion. Eventually the flavour evolution of the neutrino ensemble becomes dominated 
by coherent processes rather than decoherence-inducing collisions. This observation is of 
practical importance, because the evolution equations then reduce to MSW form. If the 
dynamics satisfies the adiabatic condition, then the evolution becomes particularly simple. 
Actually, it turns out that adiabaticity indeed holds for the parameter space of Eg. (pJ]) . 
The low temperature evolution of the asymmetry is then approximately independent of the 
vacuum mixing angle, in the small vacuum mixing angle region. Staying with our example 
of aua — Vs system in isolation, it has been computed that the 'final' value of the asymmetry 
arises at the temperature |]2B| 




T/^ 0.5 ^1^ MeV. (12) 



The magnitude of the final value was calculated to be [26|, 



Ll^ ~ 0.29/1 for |5m^J/eV2 ~ 1000, 
Ll^ ~ 0.23/1 for 3 ~ |(5m^J/eV^ ~ 1000, 

Ll^ ~ 0.35/1 for 10"^ ~ \5ml^\/eY^ ^ 3, (13) 
where h = {T^/T^Y . Similar results also hold for ordinary - mirror neutrino oscillations. 



D. The Big Bang Nucleosynthesis Epoch 



At temperatures of a few MeV, weak interaction rates start to become smaller than the 
expansion rate of the universe. This causes the ordinary neutrinos to fall out of kinetic 
and chemical equilibrium with the background plasma. It also signals the onset of the 
BBN epoch because of the end of nuclear statistical equilibrium. For the typical parameter 
space of interest in the EPM, we will show that a significant electron neutrino asymmetry is 
generated by and during during this epoch. This will have important implications for BBN, 
and one of the major goals of this paper to compute this effect. 



III. QUANTUM KINETIC EQUATIONS FOR ORDINARY-MIRROR NEUTRINO 

OSCILLATIONS 



Before we begin in earnest, we need to say something about the thermodynamics of 
the mirror particles. Because the mirror particles interact amongst themselves just like 
the ordinary particles, the mirror particles can be described by a temperature T' (and 
chemical potentials, which we assume are initially negligible). In fact, the ordinary and 
mirror particles form two weakly coupled thermodynamic systems. As in our previous 
paper [Q, we will suppose that there is an asymmetry between the temperature of the 
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mirror plasma and the temperature of the ordinary plasma so that T' <^ T. Of course if 
T' = T then a neutrino asymmetry would not be expected to develop. The energy density 
of the mirror sector would then double the expansion rate of the universe. In this case the 
reasonably successful BBN predictions would be lost. However one should remember that 
exact microscopic symmetry does not imply exact macroscopic symmetry. In reality if the 
ordinary and mirror particles are only in very weak thermal contact there is no compelling 
reason for T' = T. Note that the assumption that T' <^T does not imply that the amount of 
mirror baryonic matter in the universe today is less than the ordinary baryonic matter. The 
origin of baryon number (and mirror baryon number) is not understood at the moment, so 
no definite conclusions can be drawn regarding the amount of mirror baryonic matter (and 
hence mirror stars and so on) in the universe today. Actually there are strong astrophysical 
arguments for the existence of a large amount of dark matter in the universe, and this 
suggests that the mirror baryon number is comparable or even greater than the ordinary 
baryon number. However, as far as the earlyQ universe is concerned, the precise value of the 
mirror baryon number should be unimportant since the energy density will be dominated by 
the relativistic degrees of freedom (neutrinos, electrons/positrons and photons). Thus, when 
we use the term 'mirror matter' below, we will be referring to the 'light' mirror particles, 
that is the mirror electrons/positrons, mirror photons and mirror neutrinos, since these are 
the mirror particles which affect the expansion rate of the early universe. 

We now discuss the Quantum Kinetic Equations (QKEs) for a two-flavour subsystem 
consisting of an ordinary neutrino z/^ and a mirror neutrino i/^. We will not, in this work, 
provide an exhaustive discussion of the derivation of the QKEs or their meaning, since this 
territory is well covered in previous papers p9| , |38| -pO|] . We will, however, provide a complete 
discussion of the special features mirror neutrinos bring to the QKEs (by contrast to strictly 
sterile neutrinos). Note that two- flavour subsystems will be used as building blocks for the 
full six-flavour system in a later section. 

We will focus on evolution during the temperature regime rUe ~ T ~ m^. The plasma 
therefore consists of (i) the relativistic ordinary particles z/g, Fg, z/^, I7^, u^-, e~, and 
7; (ii) the nonrelativistic ordinary protons and neutrons (and the nonrelativistic mirror 
protons and neutrons discussed above); and (iii) whatever amount of mirror matter gets 
created through ordinary-mirror neutrino oscillations. The character of the mirror matter 
in the plasma depends on how much of it is created through oscillations. If a sufficiently tiny 
amount is created, then the mirror electromagnetic and mirror weak interactions amongst 
the mirror neutrinos will take place at a rate that is smaller than the expansion rate of 
the universe. In this case, the mirror neutrino distributions will not be of Fermi-Dirac 
form, and mirror electrons, positrons and photons will not be created. When the amount of 
mirror matter exceeds a certain level, mirror electromagnetic and mirror weak interactions 
amongst the mirror neutrinos become larger than the expansion rate. In this case, the mirror 
neutrinos produced through oscillations quickly assume a distribution of Fermi-Dirac form, 
and equilibrium distributions of mirror electrons, mirror positrons and mirror photons get 
excited in the plasma. The full plasma thus consists of two weakly coupled thermodynamic 



By 'early' we mean the time during and earlier than the BBN epoch 
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systems: the aforementioned ordinary particles at temperature T, and the corresponding 
mirror particles at a smaller temperature T'. For the case where mirror species contribute 



negligibly to the expansion rate of the universe, we have earlier shown [B^ that the inequality 



2 



must be obeyed to ensure that the mirror self-interactions are sufficiently fast to thermally 
equilibrate the mirror species. 

Our notation/convention for ordinary- mirror neutrino two state mixing is as follows. The 
weak eigenstates Ua and are linear combinations of two mass eigenstates Ua and Uh, 

Ua = cos 9a/3'l^a + 9a/3' l^b, ^'/3 = ~ siu 6'q,/j/ Z/^ + COs9a/3'l^b, (15) 

where Oap' is the vacuum mixing angle. We define Oap' so that cos29aj3' > and we adopt 
the convention that Sut^^p, = ml — m^. 

Recall that the a-type neutrino asymmetry is defined by 

^ (16) 
We also need to define an a— type mirror neutrino asymmetry, 

In the above equation, is the number density of ordinary photons. 

Note that when we refer to "neutrinos", sometimes we will mean neutrinos and/or an- 
tineutrinos and/or mirror neutrinos and/or mirror anti-neutrinos. We hope the correct 
meaning will be clear from context. 

The evolution of the ensemble of z/„ and i/^ neutrinos is described by a density matrix p„^/ 
which obeys the QKEs. A similar density matrix p^^, describes the antineutrinos. These 
density matrices pSHO] are conveniently parameterised by 



PcAp) = \[Po{p)I + P{p) ■ cr], p^p,{p) = ^[Po{p)I + Pip) ■ a], (18) 

where / is the 2x2 identity matrix, the "polarisation vector" P{p) = Px{p)^+Py{p)'y+Pz{p)2 
and a = cr^x -t- ay-y + cx^z, with at being the Pauli matrices.^ 

The quantity p is the magnitude of the neutrino 3-momentum or energy. It will be 
understood that the density matrices and the quantities Pi{p) also depend on time t or, 
equivalently, temperature T. (For the situation of negligible mirror energy density, the 
time-temperature relation for mg ~ T ~ is dt/dT ~ —Mp/5.5T^, where Mp ~ 1.22 x 
10^^ MeV is the Planck mass) . 



^Note that our previous papers used a different definition of P through the equation p = ^Pq{p)[I+ 
Pip) • a] rather than Eg. (p^) . The difference is just a matter of convention. 
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We will normalise the density matrices so that the momentum distributions of i^aip) and 
z/^(p) are given by 

NuAp) = l[Poip) + Pzip)W{p,T,0), iV,, (p) = i[Po(p) -P.(p)]iV^^(p,T,0), (19) 
where 



iV«^(p,T,/i) 



p" 



27r2 1 + exp (2^) 



(20) 



is the Fermi-Dirac distribution with chemical potential n and temperature T. Note that Pq 
is related to the total number of u^s and z/^'s of momentum p, 



N.M + K'^{P) 

iv-q(p,r,o) ' 



(21) 



while Pz{p) is related to difference, 



^ KM-K',{p) 



(22) 



Similar expressions pertain to antineutrinos. The "transverse" components Px,y{p) and 
Px,y{p) measure the degree of quantal coherence in the ensemble. Note that in subsequent 
expressions we will suppress the independent variables for notational clarity unless there is 
a chance of confusion. 

The time evolution of Pq aiid P is governed by three effects: coherent z/q, <-> i/^ oscillations, 
decoherence inducing collisions, and repopulation of z/^ and z/^ states from the background 
plasma. These effects are incorporated in the Quantum Kinetic Equations ||40|j29[| . 



dP _ dP 

'dt ~ ~dt 

dP^^dPo 
dt dt 



+ 



+ 



dP 

'dt 

dPo 



+ 



coll 



dP 

'dt 



dt 



+ 



repop 

dPo 



coll 



dt 



repop 



(23) 



where 



dP 

'dt 

dP 
'dt 

dP 

'dt 



Vap' X P, 



-DPt where Pt = -P^x + Pyy, 



coll 



(24) 



repop 



and 
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dt 







dt 



0, 



= Ru^ + R.'^. (25) 

repop 



We will explicitly define the new terms appearing above shortly. But before doing so, we 
remark that the general form of the above equations is reasonably easy to understand. 
The V X P term leads to the precession of the polarisation vector without change in its 
length. The — DPx decoherence term causes and Py to decrease in length [D > 0), 
which quantifies the rate of loss of quantal coherence. The function Ri{p) is related to 
the repopulation rate for a particle of species i with momentum p. The functions Px^y are 
unaffected by repopulation because they measure quantal coherence only. On the other 
hand Pz, being proportional to the difference in the momentum distributions of the two 
neutrino flavours as per Eg . (^2]) , receives a contribution proportional to the difference in 
the repopulation rates. The function Pq obviously remains unchanged under z/q, ^ z/^ 
oscillations, and it plays no role in quantifying loss of coherence. Since it is related to the 
sum of momentum distributions as per Eq. (pTf) , its time derivative from repopulation is 
related to the sum of the repopulation rates. 

Similar equations are satisfied for the antineutrino functions Pq and P, with the substi- 
tutions 

V„;3'-^V,^,, D^D, R^^Rj. (26) 

We now explicitly define the terms appearing in these equations. 

The function V^^/, which is related to the effective matter potential, drives the coherent 
aspect of the evolution of the density matrix. Importantly, Vq,^/ depends on the neutrino 
and mirror neutrino asymmetries. It is given by IST 



where j3 and A are 



V„^, =/3x+Az, (27) 



m = sin2^,^,, \{p) = ^[cos20„^, - h{p) ± a{p)l (28) 

in which the +(— ) sign corresponds to neutrino (antineutrino) oscillations. The dimension- 
less variables a{p) and h{p) contain the matter effects [^, being the matter potential divided 
by Sml^p,/2p. For z/q, ^ v'p oscillations a(p) and h{p) are given by 

4C(3)v/2GpT3L("/^')p 4C(3)v/2GprM.p^ 
«W = ' = 2. 2 /i^2 ' (29) 

where ({3) ^ 1.202 is the Riemann zeta function of 3, Gp is the Fermi constant, Mw is the 
VT— boson mass, Ae — 17 and Af^^r — 4.9 (for rrie ~ T ~ m^). The expression for b{p) is 
valid provided that the plasma has a negligible component of mirror energy density. The 
quantity is given by 

Li-P') = L^o.) (30) 
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where 



+ L,^ + + L,^ + r], L'(^) 



L,,i 



V 



(31) 



Recall that the term r] is due to the asymmetry of the electrons and nucleons and is expected 
to be very small, 77 ~ 5 x 10~^°. The mirror analogue, 1]' , will also be taken to be very 
small. For antineutrinos, the corresponding function Va/s' is obtained by the substitution 
) ^ — L*^"^ \ The MSW resonance conditions are given by 



0, 



(32) 



where Pres is the resonance momentum. 

The term D{p) is the decoherence or damping function. When the number density of 
mirror species is much less than the number density of ordinary species, it is given by fS^ 



Dip) 



Tip) 



(33) 



where Tip) is the total collision rate of a z/^ neutrino of momentum p with the background 
plasma.|3 From Refs. 12^,^ it is given by 



rip) = VaGlT' 




(34) 



where (p) ^ 3.15T is the average momentum of the ordinary neutrinos, i/e — 4.0 and 



l/^T- ~ 2.9 (for the rUe ~ T ~ epoch we are considering). The total collision rate for a 



v'g mirror neutrino of momentum p is roughly. 



r'(p) 



(f)'r(p) 







if T' obeys Eq.(14) 
otherwise 



(35) 



In the presence of neutrino asymmetries, the collision rates for neutrinos and antineutrinos 
differ. The collision rates quoted above hold when the asymmetries are small, with the 
antineutrino rate being approximately equal to the neutrino rate in that limit. Note that 
in the parameter space regime we are considering, neutrino asymmetries do not become 
large until temperatures are sufficiently low that collisions can be approximately neglected. 
Therefore, the dependence of the collision rates on the neutrino asymmetries is never of 
practical importance. 

The repopulation functions R^^ and R^,' are given by 



R,, 



Ru' - 



-\iP. 



P. 



2(^0 -n) 



(36) 



^'^If the number density of mirror species is significant, Dip) must also include the collision rate 
of u'pS with the background mirror particles. 
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where 



KuAP) ■ 



with /ij being the chemical potential for species i. For antineutrinos, fiu^ is replaced by fi-p^ 
in the above equation. The approximate equality sign in Eq.(p6D indicates that the right- 
hand side is not an exact result. It holds when all species are in thermal equilibrium apart 
from z/q and i^^, which are instead approximately in equilibrium. [See Ref. for the exact 
form of Eq. (PSj).] The two terms R^^ and R^'^ are due to the repopulation of z/q, states by 
ordinary weak interactions, and the repopulation of u'^ states by mirror weak interactions, 
respectively. 

In order to integrate Eqs. (|2^) and (^), we need to relate the chemical potentials ap- 
pearing in Eqs. (|36|) and ( |37D to the asymmetries appearing in Eq.(0). In general, for a 
distribution in thermal equilibrium, 

1 x^dx 1 r°° x'^dx 



4C(3) Jo 1 + e^-^^" 4C(3) Jo 1 + e^"/^^ 
where fla = f^u^/T and fia = fi-Pa/T. Expanding out the above equation 



L ~ 

24C(3) 



Tr'if^a - /i^) + Qii^i - K) In 2 + {fit- fii) . (39) 



This is an exact equation for jX^ = —fia, otherwise it holds to a good approximation provided 
that /ia a ~ 1. For T ~ T^cc^ where T^^^, ~ 2.5 MeV and T^^'^T ~ 3.5 MeV are the chemical 
decoupling temperatures, /ij,^ ~ "A^f^ because inelastic processes such as I'a^a ^ e+e~ 
and e+e^ ^ 77 are rapid enough to make fia + P^a — i^e+ + /^e- — 0. However, for 
1 MeV 

r-u ^ r-u T^^^^ weak interactions are rapid enough to approximately thermalise the 
neutrino momentum distributions, but not rapid enough to keep the neutrinos in chemical 
equilibrium.!^ In this case, the value of jla is approximately frozen at T ~ TJ^^ (taking for 
definiteness L^^ > 0), while the anti- neutrino chemical potential jla continues evolving until 
T ~ 1 MeV. For T ~ 1 MeV, the exact form for the right-hand side of Eq.(p5D should be 
used. 

The neutrino asymmetries that appear in A, R and their antineutrino analogues are 
in principle calculated from the density matrices. Recall that the neutrino asymmetry is 
defined in Eq.(|l6l). The number density of is 



.^=/ KJp= -iPo + P.)N^^{p,T,0)dp, (40) 
JO Jo 2 
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The chemical and thermal decoupling temperatures are so different because the inelastic collision 
rates are much less than the elastic collision rates. See, for example, Ref. |43| for a list of the collision 
rates. 
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so that 



2n-^ Jo 



(Po + n) - (Po + P.) N'''{p,T,0)dp. 



(41) 



Although it is in a strict technical sense redundant, it is useful to derive an equation for 
the rate of change of lepton number. It is given by 



Thus, using Eq. ([T9|) , 



dt 



2n^ Jo 



dLj^ 
dt 



dPo m 
dt 




dPo dP. 



dt 



dt 



dt 



N'^'ip, T, 0)dp. 



(42) 



(43) 



This equation can be further simplified using the QKEs and the fact that the repopulation 
does not directly affect the lepton number to obtain. 



dL^ 
dt 



2n-y Jo 



' (3{Py-Py) N'^{p,T,0)dp. 



(44) 



In our numerical work, this equation is the one actually used to calculate the lepton number 
that appears in the QKEs. Note for future reference that a limiting case of these equations 
will take centre stage when we come to study the Low Temperature Epoch. 

The last piece of information needed is the evolution equation for the mirror sector 
temperature T' . This is obtained by using a conservation of energy argument that was first 

oscillations with the mirror 



presented in Ref. [^. It goes as follows: Consider z/q, v'^ 
interactions felt by v'^ artificially switched off. The energy density of the v'^ and v'^ states is 
then given by 







2 Jo 



(N,>^+N^^^)pdp 

P^-P, + P^- p\ pN'^^ip, T, 0)dp. 



(45) 



Now switch on the mirror self-interactions. They will quickly distribute this energy density 
amongst all of the relevant mirror species: the three mirror neutrinos and antineutrinos, the 
mirror electrons and positrons, and the mirror photon. However, the energy density that is 
being fed into the mirror sector by ordinary- mirror oscillations is still given by the righthand- 
side of Eq.(^). The rate at which energy density is being transferred from the ordinary 
to the mirror sector is therefore equal to the time rate of change of the righthand-side of 
Eq.(H5]) due to oscillations only. Therefore we conclude that 



dp' 










'dt 




2 


lo di 












where p' = 3p^'^ 


+pe'+P^ 


is 


the total 



P, + Po-P^pN'''{p,T, 0)dp 



(46) 



equation for T' is obtained by combining Eq.(|^) with the cosmological red-shifting of T'. 
To this end, consider the quantity 7p where 
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P \T) 

with p = ^^T'^ being the total energy density due to ordinary species. This ratio of 
energy densities does not red-shift. Its total rate of change can therefore be calculated from 
Eqs.dll), dH) and (H) to yield 

which is the required evolution equation for T'. Note that we implicitly assumed in the above 
derivation that the redistribution of mirror energy density from i/^'s to the other mirror 
species did not affect the expansion rate of the universe. This is a good approximation 
provided that the mirror energy density is small, that is 7p <^ 1. This is generally expected 
to be the case since 6N^^eS ~ 0.6 ^ 7p ~ 0.1. 

We end this section with an example of the evolution of neutrino asymmetry generated 
by two flavour ordinary - mirror neutrino oscillations. In Fig.l the evolution of the r- 
like asymmetry is plotted. For our present illustrative purpose, we have considered the 
evolution of L^^ under the influence of the Ur — i''^ oscillation mode only. The parameter point 
Sm^^i = —50 eV^ and sin^ 26r^' = 10~® has been chosen. The initial L^^ is set to zero. Notice 
that L^^ evolves from zero to a value which approximately cancels the baryon asymmetry 
by T ~ 70 MeV. The asymmetry then remains constant until the critical temperature 
Tc ~ 38 MeV when explosive growth begins. Shortly thereafter, the explosive growth phase 
gives way to power law growth. During the power law phase, the High Temperature 
Epoch evolves into the Low Temperature Epoch. Of course the behaviour shown in Fig.l 
is quite general and in fact quite similar to ordinary - sterile neutrino oscillations. The 



latter have already been studied in some detail in previous papers 1 24 -26, 28, 29 1. Finally 



note that we have plotted \Lty^\. This is because L^_^ changes sign at the critical temperature 
(the reason for this behaviour has been discussed in Ref. |2^). For values of sin^ 2^ large 
enough, our numerical results indicate that the sign of L^,^ initially oscillates and thus the 
final sign of the asymmetry may be random. This may lead to different regions of space 
having different neutrino asymmetries (as suggested earlier in Ref. p4|). Also it should 
be mentioned that the effect of statistical fluctuations on lepton number asymmetry is an 
important open problem, and consequently it is also possible that the sign of the asymmetry 
may turn out to be random even for small values of sin^ 29. For the purposes of the present 
paper, we acknowledge the indeterminate nature of the asymmetry by considering the two 
possible signs in all our numerical work. 

In the next section we will discuss the High Temperature Epoch in the EPM. Our main 
goal there will be to demonstrate the consistency of the z/^ —>■ u'^ solution to the atmospheric 
neutrino problem with BBN for a range of parameters. 



IV. HIGH TEMPERATURE EPOCH: CONSISTENCY OF u'^ SOLUTION OF 

THE ATMOSPHERIC NEUTRINO ANOMALY WITH BBN 

Previous work has shown that the QKEs for ordinary-mirror (or ordinary-sterile) oscil- 
lations will imply the explosive creation of neutrino asymmetries provided some fairly mild 
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restrictions on the oscillation parameter space are imposed. In particular, the involved 
must be negative, and the vacuum mixing angle Oap' must be in the approximate range 

10-^° ^ sin^ 2e^p, ^ few x 10"^ . (49) 

The lower bound comes from the requirement that the oscillation mode be sufficiently 
strong]^, while the upper bound was derived in Ref. from the requirement that Va ^'13 
oscillations, considered in isolation, not spoil BBN [Eq.(^) takes for definiteness that 
SN,y^eg ~ 0.6]. Once created, the large neutrino asymmetry or asymmetries will suppress 
other ordinary-sterile oscillation modes for a range of parameters. 

For the generic parameter space region considered here (see Sec. II), the oscillation modes 

^'f,, i^'e, -* (50) 

could all satisfy the above criteria. We will call these the "lepton number creating modes" . 
The other oscillation modes, including the z/^ — >■ u'^ mode that hypothetically solves the 
atmospheric neutrino problem, tend to destroy a linear combination of asymmetries. 
Lepton number amplification begins at a critical temperature T^, given roughly by 

^.«l<-^)'MeV. (51) 

where the oscillation parameters pertain to the two-fiavour lepton number creating mode 
responsible. The mode with the largest |5m^| will therefore be expected to create lepton 



number first, provided its vacuum mixing angle is in the range of Eq.(|U). Within the 
scenario of Sec. II, the Uj. — >■ u'^ and — >■ z/g modes will have the largest 6m'^ values. Which, 
if either, of them dominates lepton number creation then depends on their specific oscillation 
parameters. For the sake of a plausible example, we will suppose that Ur v'^ dominates, 
even though Vr ^'e has a slightly larger brv? . This is because we expect Qre' ^ Qrii' as per 
Sec. II. Basically, we will work in the parameter space region where Ore' is negligibly tiny. 
So, we are led to consider the four fiavour subsystem. 



^ ^'^ 



: X : ■ (52) 

^ V'^ 

Further, we decompose this four fiavour system into the two fiavour subsystems indicated 
by the arrows above. Some discussion of the justification for this sort of decomposition 



can be found in Ref. p9|. Heuristically, it is expected that this simplifying assumption is 
justified because the MSW resonance momenta of each of the oscillation modes are generally 
different. The ordinary-ordinary and mirror-mirror modes are governed by the same mixing 
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But note that the vacuum oscillation amplitude can still be tiny! 
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angle ^^t, the modes z/^ z/^ and z/^ u'^ are both governed by another mixing angle ^t^/, 
whilp tVip 7y_ .. ^ 7/^ modes are maximally mixed. The squared mass difference ^m^^/ is 



set by the atmospheric neutrino data to be in the range quoted in Eq.(|^. The other mass 
parameters, 5m1^, and 5m^^,, are free, subject to the restrictions discussed in Sec. II and 
summarised in Eq.(|^. 

For simplicity and the sake of the example, we will set |(5m^^/| to be so small that the 
associated oscillation mode can be neglected.Q The precise value of the mixing angle O^t 
is unimportant, provided it is small. The <-> z/^ mode has almost no effect until lepton 
number is large (here 'large' means greater than about 10~^), because of the approximately 
equal number densities of the two species involved. In the High Temperature Epoch being 
considered in this Section, lepton number will always be small. For a given Sm?^^, we are 
therefore effectively left with two free parameters: 6*,-^' and Sm"^^,. Our task is to find the 
region of this parameter space for which the z/^ —>■ v'^ solution to the atmospheric neutrino 
problem is consistent with BBN. For the EPM[^this calculation was first performed in Ref. 
|3^ within the static approximation. We improve on this approach here through the use of 



the QKEs. 

We now write down the equations we must solve. We begin by introducing three two- 
flavour density matrices, 

Pr/.' = ^(i^O + O- ■ P), Pm/^' = ^(<5o + O" ■ Q), P;.r' = ^(5*0 + Cr ■ S), (53) 

for the three significant oscillation modes, 

T^r ^ l^'f,, i^t, ^ ^ K, (54) 

respectively. Since u'^ is common to the first pair of modes, and z/^ is common to the second 
pair, we have the constraints 

N^' 1 1 

-(Po-p.) = -(go-g.), 



A^«=i(p,T,0) 2' ' 2 

N 1 1 

7:iQo + Qz) = -iSo + S,). (55) 



iV<'q(p,T,0) 2'^ ^ ' 2 

Extending the two flavour case discussed in the previous section, the time derivatives of 
the functions Pq, Qq, Sq, P, Q and S are observed to receive contributions from each of 
the three oscillation modes, from decohering collisions, and from repopulation. Denoting a 
generic function by F, we have that 



^^This will be the case provided that |(5m^^, | ~ l^^fif^'l- If I'^'^tt'I much larger than \Sm'^^,\ 
then the resulting 'allowed region' will be significantly reduced. 

For the case of strictly sterile neutrinos, this calculation was done in the static approximation 
in Ref. [^] and by numerically integrating the quantum kinetic equations in Ref. [p^ . 
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OF _ OF 
'dt~~dt 



+ 



OF 



dt 



+ 



dF 



dt 



+ 



dF 



dt 



+ 



dF 



coll 



dt 



repop 



From the two-flavour formalism described in Sec. Ill we have that 

= V.^, X P - DPt, 
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(9P 


'dt 




'dt 


dQ 
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dQ 


dt 




dt 






dS 




+ 


dS 


'dt 




'dt 



coll 



coll 



coll 



X Q - DQt, 



\.r' X S - DSt. 



It is also clear that 



dPo 




dPo 




dt 




dt 


coll 


dQo 




dQo 




dt 




dt 


coll 


dSo 




_ dSo 




dt 




dt 


coll 



= 0, 
0. 



We also obviously know that 



dP 


dPo 




'dt 


dt 




dS 


dSo 




'dt 


dt 





0, 
0. 



(56) 



(57) 



(58) 



(59) 



Consider now the contribution of z/^ v'^ oscillations to the evolution of P and Pq- First 



of all, the transverse components P^ y receive no contribution, 



dP, 



x,y 



dt 



(60) 



because they are affected only by decohering collisions. The evolution of Pz and Pq can be 
obtained from Eq.(^) by noting that 



d_ 
di 



{Po - P. 



d_ 

dt 



(Qo - Q. 



dQz 



dt 



(61) 



SO that 
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dPo 
dt 

dt 



+ 



2 dt 
2 dt 



(62) 



having used Eq.(0). The expression for {dQz/dt)\^ is obtained from Eg. {^7\) . Finally, 
we have that 



dt 



0, 



rcpop 



ropop 



repop 



r' 



(63) 



This completes the specification of the evolution equations for P and Pq. 

The evolution of S and 5*0 due to — > u'^ oscillations and repopulation is handled in a 
very similar manner to yield, 



dS. 



dt 
dSo 



dt 

dS^_ 



dt 

dSx^y 



dt 



repop 



repop 



repop 



0, 



^1 dQ, 



2 dt 
dQ, 



2 dt 



0, 

r 
r' 



Ky'^ - -(5*0 - Sz) 



(64) 



The completes the specification of the evolution equations for S and Sq. 

Finally, we have to specify the time rate of change of Q and Qq under the influence of 
Vr ^ ^'^ oscillations, z/^ ^ u'^ oscillations and repopulation. We note first that Qx,?/ are 
unaffected by these processes: 



dQ 



x,y 



dt 



dQx,y 



dt 



dQx,: 



dt 



0. 



(65) 



rcpop 



Then, from Eq.(ES|) we see that 
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dQo 



dt 



dt 
dQo 



1 d_ 

2 di 

1 d 



{Po - Pz 
{Po - P. 



1 dP, 



dt 



2 dt 



2 9t 
^1 



2 dt 
1 9S 



2 9t 

1 dS-. 
+7: 



Finally, 



d_l 

dt2 

d_l 

dt2 



{Qo + Q.) 
(Qo - Q.) 



rcpop 



repop 



r 
r' 



2 dt 



K,--{Qo + Q.) 



(66) 



(67) 



specifies the remaining repopulation equations. 

This completes the list of quantum kinetic equations for our system. Of course, because 
of the constraints in Eq.(^), some of these equations are redundant. 

To use these equations one needs, (i) equations connecting neutrino asymmetries with 
chemical potentials, and (ii) an evolution equation for the temperature T' of the mirror 
plasma. The chemical potentials are calculated in exactly the same way as discussed in 
Sec. III. The T' evolution equation is obtained by extending the two-flavour derivation of 
Sec. Ill in the obvious way. As before, it is best to first imagine that the mirror electroweak 
interactions are artificially switched off. The energy density in mirror states is then entirely 
due to the u'^ and u!^ species: 





1 

2 Jo 



N,,^ + Nt,>^ + N,, + Nt,^^) pdp 
Po-P, + Po-Pz + So~S, + So- S^ pN'^'ip, T, 0)dp. 



(68) 



With mirror electroweak interactions now switched on, this energy density is distributed 
amongst all of the relevant mirror species. The rate at which energy density is being trans- 
ferred from the ordinary to the mirror sector is thus 



dp' 
It 



d 



P,-P, + P,-P^ + So-S, + So- 5 J pN'\p, T, 0)dp, (69) 



2 Jo dt 

where p' is the total energy density in mirror species, and 



d_ 

dt 



d 
dt 



+ 



d_ 

di 



+ 



d 
di 



(70) 



Introducing •jp as per Eq.(^7D and using the QKEs we obtain, 
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dip ^ l_ 

dt 2p 



r [Pr^.' {Py + Py) + {Qv + Qy) + ^^r' [Sy + pN^'^p, T, 0)cip (71) 

'J 



as the T' evolution equation, where 

= ^^sin20.^,, /?^^, = ^^, = _^sin20.^,. (72) 

Note that ~ —i^iir' because 5m^^, ~ —6m'j^^, for the parameter space of interest. 

We first present the main result of numerically solving the above equations. After doing 
so, we will provide a physical description of what lies behind the mathematics. The main 
result is displayed in Fig.2, which shows the region of (5m^^/, sin^26'T-^/) parameter space 
which is consistent with Big Bang Nucleosynthesis, for various values of 5m^^; motivated 
by the atmospheric neutrino anomaly. The allowed region lies above the relevant solid line 
(which corresponds to a particular 5rn?^^,), and to the left of the dash-dotted line. The solid 
lines arise from solving the QKEs, while the dash-dotted line is the upper bound quoted in 
Eq.(^9D applied to the lepton number creating mode z/,- ^ z/^. For the sake of definiteness, 

we have adopted 5Ny^cE ~ 0.6 as the BBN bound on the expansion rate of the universe 
(expressed as an equivalence to additional relativistic neutrino flavours, as is customary). 
Of course, at the present time there is some confusion regarding the value of this bound, 
due to conflicting primordial element abundance measurements. The value of 0.6 was chosen 
for illustrative purposes only. The position of the dash-dotted line depends on the SNy^^s 
chosen. The solid lines, on the other hand, define sharp transition regions. Below the lines, 
the mirror sector comes into thermal equilibrium because of the eventual copious production 
of v'^ from i/^ — > v'^ oscillations. Above the lines, essentially no mirror matter is produced by 
this oscillation mode. Of course, the closer one gets to the dash-dotted line, the more mirror 
matter is produced by the Vr — ^ v'^ mode. For sin^ 20r^i ~ 10~^ the results obtained here 
using the QKEs are almost identical to those obtained earlier using the static approximation. 
The results differ at large values of the mixing angle mainly because lepton number is created 
rapidly enough to spoil the validity of the static approximation.]^ 

The importance of Fig.2 lies in its demonstration that the —>■ v'^ solution in the EPM to 
the atmospheric neutrino problem is cosmologically consistent for a large region of oscillation 
parameter space. Furthermore, most of this region sees the u^- having a cosmologically 
interesting mass. In particular, note that the hot dark matter region marked as a shaded 
band in Fig.2 has a significant overlap with the allowed region from BBN (we will discuss 
more about the dark matter region in section VI). It is also interesting to note that the 
BBN allowed region implied by the z/^ v' solution to the atmospheric neutrino problem 



is larger than the corresponding region p5| , p8| obtained when v'^ is replaced by a strictly 



sterile neutrino. The production of sterile neutrinos tends to delay the onset of the rapid 
exponential growth |^ which means that by the time it occurs the z/^ ^ z/^ oscillations 



^^A complete discussion of the static approximation can be found in Refs. |25,29|. In particular, it 
was shown in Ref. ||2^ that the static approximation is an adiabatic-like approximation for partially 
incoherent oscillations in the small vacuum mixing angle parameter space regime. 
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can destroy L^f^f^"^ more efficiently since they are not damped so much by the colhsions. 
In the case of the mirror neutrino z/' the mirror electroweak interactions have the effect 
of reducing their number density, so the production of lepton number is not delayed. Of 
course a Ur in the eV mass region is currently being searched for in the short baseline 
Nomad/Chorus experiments. Such experiments are extremely important to test for the eV 
tau neutrino which is suggested by Figure 2. Unfortunately, we cannot predict sin^ 29^-^, so 
these experiments will either discover z/^ ^ z/^ oscillations or constrain sin^ 29r^- 

Before closing this section, we will discuss some of the numerical details of performing 
the above computation, which will entail also a discussion of the physics of the result. 

MSW resonances play a key role in the evolution of the system. It is instructive to 
examine the connection between the neutrino asymmetries and the resonance momenta 
of the three important two-flavour modes within our system. Note, first of all, that the 
effective potentials of the three modes depend on different linear combinations of neutrino 
asymmetries: 

L(^'^') = 2L,^ + 3L,^-L,,, 
L('^'^') = 3L,^ + 2L,,-L,,, 

L(^^')=4L,^ + 3V + L,,, (73) 
Lj,^ = 0, and we have used conservation of lepton number^ 



where we have set L, 



Lu^ + + Ly^ + L^i^ + L^i + Lyi^ — 0, 



(74) 



to eliminate L^i . The lepton number creating mode Vr ^'n generates a nonzero L,j^ , which 



means that L is also nonzero. The latter quantity then suppresses 



oscillations. 



provided that it grows sufficiently quickly for a sufficiently long period of time. This is not 
inevitable, because the effect of the lepton number destroying mode z/^ — » u'^ is to try to 
destroy L'-'^^'-' through the creation of nonzero values for Ly^ and Lyi^ to compensate the 
nonzero Ly^ and L^/^ produced by the Vr v'^ oscillations. The essence of the calculation 
presented above is the determination of when L'^^^^^'^ is driven to zero, and when it is not. 
This depends on the oscillation parameters, as summarised in Fig.2. 

The resonance momentum for the Vr ^'^ mode and its antimatter analogue are given 

by 



T 



T 



6or2 boT^ 



/aoL(V)y ^ 4|(5m2^,|cos2^^^, 



+ 



aoL(^'^')y 4|5m2 ,|cos2^^^ 



+ 



(75) 



Of course the sum of lepton numbers need not be exactly zero. However, we can set the sum 
to zero without loss of generality provided that the sum is not large. 
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The p- and T-independent quantities ao and 60 are defined through, 



= -''o^Zr2 ' Kp) = -boT-^, (76) 



T3L(^f^')p TV 
a[p) = —a^ ^ — 

leading to 



ao = 5 , &o = TTT2 ■ (77) 

In the following we will consider the L'^'^^'^ > case for definitenessQ Before the rapid 
exponential creation of lepton number (that is for T > T^), the neutrino and antineutrino 
resonance momenta for the lepton number creating modes are equal. As L^'^^ •* gets expo- 
nentially created, the neutrino resonance momentum p,-^/ moves rapidly to infinity, while 
the antineutrino resonance momentum p^^, remains at a value of order T. Numerical calcu- 
lations show that p^^f /T typically takes a value in the range 0.2 — 0.6 at T ~ Tc/2 where T^. 
is the critical temperature at which lepton number creation begins. These observations are 
important, because they mean that v'^ oscillations are unimportant after the creation 

of lepton number, while v'^ oscillations remain very important. This is simply because 
the neutrino resonance momentum has moved to the tail of the Fermi-Dirac distribution, 
while the antineutrino resonance momentum is within the body of the distribution. (Note 
that these observations will play a central role in the next section.) 

The resonance momenta for the z/^ — > z/^ and — > V'^ modes are given by 

P^^a_^ ^ = ifL(-')>0; 
T boT^ T 



if L^'^'^') < 0. (78) 



In the region of parameter space where L'^^^^^"> is not driven to zero, we see that p^^i/T gets 
driven to infinity (staying with the Lj,^ > case), while p^^, /T stays at zero. In the region of 
parameter space where L^'^'^'^ gets destroyed, we see that p^i^j.' /T moves from zero to a finite 
value as L^^ gets created, then moves back towards zero as the compensating L^^ is induced. 
There is a sharp transition between these two possibilities for the evolution of p^^^/fT, with 
the boundary given by the solid lines in Fig.2. Above the solid line L^'^^ •* is created early 
enough and is large enough so that the z/^ is never significantly populated by z/^ z/^ 
oscillations, these oscillations being heavily suppressed by the matter effects resulting from 
the large L^^f^ \ Below the solid line the the z/^ states would eventually become populated 

by z/^ <-> z/^ oscillations in the temperature range 6 ~ T/MeV ~ 10. Furthermore, the other 
mirror particles would also become populated due to the mirror weak interactions, which 
would effectively double the energy density of the universe prior to the BBN epoch. 
The z/^ — >■ ly'j. and z7^ — > V'^ resonance momenta are given by 



^"^A discussion of the overall sign of the asymmetries created can be found in Ref. [25|. 
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T 



P. 



fJ,T' 



ao 



4|5m2 ,| cos 26*, 



± 



6oT2 \ 1^ boT^ 



A\6ml^, \ cos2e. 



(79) 



Because the sign of is positive, we see a qualitatively different behaviour for these 

resonance momenta compared to their mirror reflections in the i>r + z/|j subsystem. Before 
the creation of lepton number, there are no solutions to the resonance conditions. If lepton 
number evolves to the point where 



ao 



boT^ 



baT<^ 



(80) 



then (taking the L^^'^''^ > case) z/^ ^ v'^ comes on resonance at a finite value of the 
momentum, while ^ V'^ never comes on resonance. In the region of parameter space 
where L^^ dominates, it is easy to show that this point occurs when 



T 



3 T 



2. 



(81) 



Thus by the time sufficient Ly^ and L^i^ asymmetries have been generated by 
oscillations for the v'^ oscillations to have a resonance momentum, the i^'^ - - -^ii 

oscillation resonance momentum is already into the tail of the distribution (for the parameter 
space where negligible L^^ has been created by z/^ ^ u'^ oscillations) .[^ For low temperatures, 

T ~ Tc/3, the 6— term in the matter potential can be approximately neglected and the 
resonance momentum can be derived from a{p) ~ cos 26'^^/, leading to 



T 



6ml^, cos 29 



(82) 



Observe that the effect of the ^ v'^ oscillations is to decrease |L*^'^^')| and hence to increase 
p^r'/T. By the time T ~ Tc/2, the resonance momentum has reached p^^r'/T ~ 15. We 
mention this here because it will be important in the following section. 

The observations about the evolution of resonance momenta made above are relevant 
to the numerical integration of the quantum kinetic equations. Because this integration 
is CPU time consuming, we employ the useful time saving approximation of integrating 



Actually, in the psuedo-Dirac case |33| where there are no mirror interactions the + ly'^ 
oscillation system is more important. The reason is that, in the psuedo-Dirac case, these oscillations 
destroy exactly the same combination of lepton numbers as does the mode, that is, in this 

case L^^'^'^ = L^f^^'\ Thus, in the pseudo-Dirac alternative to the mirror scenario, the ^ v'^ 
oscillations can help the v'^ oscillations destroy L'^^^'\ which means that the 'allowed region' 

can be significantly reduced. 



27 



the oscillation and collision driven aspects of the evolution in the region around the MSW 
resonances. Since the precise details and justification of this have been covered in Ref. [28|, 
we will not repeat the discussion here. 



V. LOW TEMPERATURE AND BBN EPOCHS: EFFECT OF OSCILLATIONS 

ON LIGHT ELEMENT ABUNDANCES 

A. Introduction 

The primordial deuterium to hydrogen (D/H) ratio can be used to give a sensitive 
determination of the baryon to photon ratio r) which, given the estimated primordial ^He 
mass fraction, can be used to infer the effective number of light neutrino flavours A'^i.^efr 
during the BBN epoch. This value can then be compared with the predictions for Ni^^eS 
from various models of particle physics to flnd out which ones are compatible with standard 
BBN. For example, the minimal standard model predicts N^^^g = 3. At the present time. 



most estimates favour iV^^eff < 3.6 and some estimates favour A^'j.efr < 3.0 |Q. Of course, 
even if a model of particle physics is shown to be incompatible with BBN, this does not 
necessarily mean that the model is incorrect, since it is also possible that one of the standard 
assumptions of BBN may not be correct . 

For gauge models with mirror or sterile neutrinos, one in general expects N^^^s 7^ 3. In 
fact, A^'j.eflf niay be less than three or greater than three. The prediction for A^^^efr depends 
on the oscillation parameters in a given model. One possible consequence of ordinary-mirror 
(or ordinary-sterile) neutrino oscillations is the excitation of mirror neutrino states, which 
typically leads to an increase in the expansion rate of the universe and thereby also increases 
Ny^eS- Another possible consequence of ordinary-mirror neutrino oscillations is the dynamical 
generation of an electron-neutrino asymmetry. This also has important implications for 
BBN, as it directly affects the reaction rates which determine the neutrino to proton {n/p) 
ratio just before nucleosynthesis. If the electron neutrino asymmetry is positive then it will 
decrease A^^,efr, while if it is negative then it will increase Ny^es- 

The neutron to nucleon ratio, is related to the primordial Helium mass fraction, 

l^p, byO 

Yp = 2X„ (83) 

just before nucleosynthesis. The evolution of X„(t) is governed by the equation, 

rlX 

_^ = _A(n p)X„ + \{p -.n){l- X„), (84) 
where the reaction rates are approximately 



X{n — i> p) ~ \{n + ^ p + e ) + X{n + ^ p + 
\{p n) X{p + e~ ^ n + Ue) + X{p + Ve ^ n + e"*"), (85) 



^^For a review of helium synthesis, see for example Ref. |46] 
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depend on the momentum distributions of the species involved. The processes in Eq. (^) for 
determining n ^ p are only important for temperatures above about 0.4 MeV. Below this 
temperature the weak interaction rates freeze out and neutron decay becomes the dominant 
factor affecting the n/p ratio. An excess of over Fg, due to the creation of a positive L^^ 
would change the rates for the processes in Eq. (pS]) . The effect of this would be to reduce 
the n/p ratio, and hence reduce Yp. Neutron decay is not significantly altered by lepton 
asymmetries. It is quite well known that a small change in Yp due to the modification of 
and T7g distributions does not impact significantly on the other light element abundances 
(see for example Ref. |^). A small modification to the expansion rate, using the convenient 
unit N^eS, primarily affects only Yp, with 



SYp ~ 0.012 X SN^^eS- (86) 

In Appendix A we describe in detail how we compute the effect on Yp due to the modified 
Ue and I/g distributions. 

In two previous papers [^,^, we studied the implications for BBN of oscillations within 
two distinct four-neutrino-flavour models which featured the three ordinary neutrinos and 



one sterile neutrino. In Ref. |2g], a model with the mass hierarchy m^^^ ^ mu^,m^^,miy^ 
was considered. In this case ^ T^s oscillations resulted in an excess of Ur over (in 
the case where L^^ > 0) thereby generating a large tau-neutrino asymmetry. It was shown 
that if Ft- — Fg oscillations also occurred, then some of the tau-neutrino asymmetry was 
reprocessed into an electron-neutrino asymmetry. The effective number of neutrino flavours 



found in Ref. ||26| was either 2.5 or 3.4, depending on the ambiguity for the sign of the 
asymmetry and hence the prediction for N^^^^f^ (Ref. ||2^ discusses the sign ambiguity issue). 
For a positive asymmetry, SN^^eS — —0.5 was obtained over a range of mass differences 
l^m^^l ~ 10 — 1000 eV^, while for a negative asymmetry the result was SNi^^es — +0.4. 
Later, in a separate paper with Bell [^, we considered another four-neutrino model where 



Ut- and were taken to be approximately maximal combinations of two nearly degenerate 
mass eigenstates, ui and 1^2, with m^^^rriy^ ^ m^^,m^^. In that case, we found N^^^s ~ 2.7 
or 3.1 depending on the sign of the asymmetry. 

As the above paragraph illustrates, the prediction for Ni,^es is a model dependent quantity. 
In the next section we will estimate N^^es in the EPM for various illustrative parameter 
ranges. 



B. Low temperature neutrino asymmetry evolution in the EPM: Case 1 

We now study the "low temperature" evolution of the number distributions and lepton 
numbers in the EPM. As discussed in Sec. II, by "low temperature" we mean the regime suc- 
ceeding the exponential growth epoch. In this regime, the evolution of the neutrino ensemble 
is dominated by coherent effects, because the decrease in the damping function D renders 
negligible the decohering effect of collisions. Repopulation, however, is still important. 

Consider, for the moment, two-flavour small angle ordinary-mirror oscillations z/q, ^ z/^. 
In the case of undamped evolution, we know from numerical integration of the exact quantum 
kinetic equations that the adiabatic approximation is valid provided that sin^ 2^q,/3' ~ 10~^°. 
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Now, coherent adiabatic MSW transitions completely convert Ua ^ v'p at the resonance mo- 
mentum of these states. For adiabatic two-flavour neutrino oscillations in the early universe 
it is then quite easy to see that the rate of change of lepton number is governed by the 
simple equation |2^ , 



'^^Vg 

dT ~ dT 

where 



-X 



d{pr:es/T) 



dT 



(87) 



X = - (N.g - N,,) , (88) 



and the case Lj,^ > has been considered (so that the resonance occurs for antineutrinos). 
Equation (|87|) relates the rate of change of lepton number to the speed of the resonance 
momentum through the neutrino distribution. Reference |^ provides a detailed discussion 
of how this equation can be derived from Eq.(^^ for the case of adiabatic evolution with a 
narrow resonance width.[] Equation (0) can be simplified using 

dT dT dU dT ' ^ ^ 

from which it follows that 

dLy^ _ dLy,^ _ fX^^S^ _ -AfXp^jT^ 



dT dT 1 fv d{Pr../T) 1 yxp; 



(90) 



T[L(<»)-L'{/3)] 



where / = 1 for d{pres/T)/dt > (that is for d{pj.cs/T)/dT < 0) and / = — 1 for 
d{pvcs/T)/dt < 0. For the multi-flavour case under analysis, coupled equations based on 
Eq.(|nD will be used. 

Of course the evolution of the lepton number can also be described using the QKEs. As 
mentioned above, they give the same answer provided that the evolution is adiabatic. In the 
case of non-adiabatic evolution, the QKEs should be used instead of the simple Eq.(^). For 
our study of the low temperature evolution of the number distributions and lepton numbers 
in the EPM, we will make use of the adiabatic approximation encoded in Eq.(0). In fact, 
it turns out that the evolution of the system in the EPM model is quite complicated. For 
example, three-flavour effects cannot be ignored, so solving the problem using the Quantum 
Kinetic Equations would be extremely complicated and (CPU) time consuming. 

We first consider the parameter region 

my.r+ - m^-r- > my^+ , ^v^- , ^v,+ , ^ve- 5 (91) 

with 



Note that when collisional decoherence in neglected, the QKEs produce standard Schrodinger- 
like MSW evolution with repopulation effects added via a Boltzmann approach. 
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< 1 eV. (92) 

We will call this "Case 1" . (Later on we will consider another case, Case 2, where m^^^ ~ eV 
as suggested by the LSND results.) In Case 1, the following oscillation modes all have 
approximately the same |5m^|, which we denote as 5mfa^j.ge- 

Vr ^ l^'e, l^T ^ l^'r ^ ^e, v'^ ^ ^^a*, 

^ Ve, Vt ^ v^,, v'^ ^ Z/g, v'^ ^ v'^. (93) 

Note that ^ivn^^^^^ ~ fnl^^. All the other oscillation modes have much smaller values. 
In fact, for Case 1, we will consider the parameter space region where the Sm^ values of all 
the other oscillation modes are small enough so that they can be approximately neglected 
for temperatures T ~ 0.4 MeV. This last condition means that these modes will not affect 
the neutron/proton ratio and hence cannot significantly affect BBN. 

In the following discussion we consider the case L,y^ > for definiteness. This means that 
the Vt- z7g and 17^ ^ v'^ oscillations generate L^^ while the other oscillations reprocess some 
of this asymmetry into other flavours. Of course a crucial issue for BBN is to find out how 
much of this asymmetry is reprocessed to the electron neutrinos, and at what temperature 
this occurs. 

In order to use Eq. (|90D , we have to employ the resonance conditions to determine the 
resonance momenta as functions of temperature and the neutrino asymmetries. We begin 
by noting that the sum of the ordinary and mirror lepton numbers are conserved by the 
oscillations, and we will suppose that they sum to zero0 

L,^ + L,^ + L,^ + L,, + L,, + L,, = 0. (94) 

Furthermore, we take as initial conditions that 

L^,^~L^^, L^^~L^,. (95) 

We will show that this assumption is robust shortly. With the above initial conditions, it 
follows that the eight oscillation modes of Eq.(^) can be classified together into four groups 
of two, each group having approximately the same resonance momentum: 



Group 1: 


Vr ^ 


— / 






Ptcs 




Group 2: 


<^ 




<^ 




Pres 




Group 3: 










Pres 


= P3, 


Group 4: 




— / 






Ptcs 





In the Low Temperature Epoch, the 6— term in the effective potential can be approximately 
neglected because it decreases as T^. This means that the resonance condition is approxi- 
mately a{p) = ±cos2^^ ~ ±1 for small angle oscillations. For ^ I'p small angle oscilla- 
tions, the resonance momentum is therefore to a good approximation given by 



Of course our results do not depend significantly on this assumption. For example, if we put 
the sum in Eq.(p^) equal to a number of the order of the baryon asymmetry then the resulting 
analysis will change very Uttle. 
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where we have used the notation defined earher in Eg. (|76D . For the four groups of oscillation 
modes in Eq.(^), 



T aoT^L, 



5 



where z = 1, . . . , 4 and 



L, = L'(-) - L'(^) = L,. - L,, = h,^ + + h,,. (99) 

Note that Eqs.(|9^) and (|95|) have been used in the above equation to express Lj,^, L^'^ and 
Li,/ in terms of the Ly^, Ly^ and Lj,/.. 

In the following discussion we will focus on the parameter space region where all of the 
oscillations are approximately adiabatic. This is extraordinarily helpful, because adiabatic 
transitions are independent of the vacuum mixing angles (so long as the mixing angles are 
much smaller than one). This means that generic outcomes can be calculated for a reason- 
ably large range of parameters, rather than having to consider small points in oscillation 
parameter space on a case by case basis. As noted earlier, two fiavour subsystems in this 
epoch of the early universe evolve adiabatically provided that the relevant sin 
This is not a very stringent requirement. In particular, there will be a large range of param- 
eters where the evolution is both adiabatic and satisfies the experimental and cosmological 
constraints. Given the parameter region of Eqs.(|9TD and (p^), consistency with BBN con- 
strains sin^ 26t-^' to be0, 



eV' 

sin^2^^^/~ few x 10"^ ( ^^^^3 — ) , (100) 



large / 



while Nomad/ Chorus constrain sin^ 26'^^ to be 



sin^ 26^^" ~ 10"^ for 6m\^^^^ ~ 40 eV^. (101) 

We now discuss the effects of each of the four groups of modes in Eq. (|96|) . 

1. The Vt- I7g and v'^ Group 1 modes: These modes have negative Srn^ values 

and thus create the relevant lepton numbers Lj,^ , L^i^ and L^i . It is important to 
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Note that sin'^ 20,-^/ = sin^ 2^t-'„ from the parity symmetry. 
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understand that if these two modes have shghtly different resonance momenta, say P" 
and Pi, then they generate lepton numbers so that Pf — ^ P^. This is tantamount 
to ensuring that the initial conditions of Eg. (pS]) hold, provided that the difference 
between the initial values of L^^ and L^^ is not too great. To see that the resonance 
momenta are dynamically driven to coincide, assume that Pf > Pi This means that 
the Vr ^ z7g resonance momentum preceeds the Vr ^ V'^ resonance momentum. Now, 
the Vr ^ V'^ oscillations convert all of the resonant V^s into z/'^'s. The closely following 
Vt V'^ resonance has a much weaker effect, since there are no 17,- 's left to convert 
into V'^ states. (For this to be true the resonance momenta must be close enough so 
that the converted Vj. states do not get completely refilled by the weak interactions 
before the trailing resonance momentum P^ passes their momentum value.) Because 
of the disparity in raw material for processing, Lj,/ is created much more rapidly then 
Lyi^. According to Eg . (p8|) , this in turn means that P" increases more slowly relative 
to P^ and thus P^ P". Obviously, if we had started with P^ > P^, then we also 
would have found that the evolution of lepton numbers is such that P" P^. Because 
the dynamics drives the two resonances in this group to approximately coincide, the 
system cannot be described in terms of two-flavour oscillations. Instead, three-flavour 
effects amongst 17^, V'^ and V'^ effect the adiabatic conversion 



(102) 



This means that as Pi sweeps through the momentum distribution, 

1 [N-,(PA AL, (PA 



1 
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1 
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1 r 
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2 2 
N-APi) + N^dPi) 



(103) 



In our numerical work the continuous momentum distribution for each flavour is re- 
placed by a finite number of 'cells' on a logarithmically spaced mesh. As the momentum 
Pi passes a cell, the number density in the cell is modified according to Eg.( |103| ) .p^ 
Of course weak interactions will repopulate these cells as they thermalise the neutrino 
momentum distributions. We will discuss this later. 

2. The v'^ ^ Vf, and v'^ ^ Group 2 modes. These modes have positive Sm^ values. 
At guite high temperatures, where the Group 1 oscillation modes are exponentially 
creating L^^, the Group 2 oscillation modes generate L,^/^, L^^ and L,,^ such that 



"Note that it is legitimate to consider probabilities, as encoded in the number density distribu- 
tions, rather than probabihty amplitudes in effecting the conversion. This is because fully adiabatic 
transitions are sufficiently 'classical'. 



33 



L'^'^'^ — L^*^^ and L'^'^^ — L^^^ 0. As already discussed in section IV, this makes 
P2/T ^ 1 a.t the onset of the low temperature epoch. Our numerical work shows that 
the initial value of P2/T is typically about 15 and decreasing by the time T ~ Tc/2. The 
subsequent evolution of P2/T is a little complicated, but can be roughly understood 
from Eqs. (|89|) and (pSf), which combine to produce 



d{P2/T) 
dT 



P2 
' T 



■4 I dL2 
T ^ T2IIT 



(104) 



or equivalently. 



diP2/T) 
dt 



P2 
T 



5.5T^ 
Mp 



1 dL2 



(105) 



By the time T ~ Tc/2, the group 2 modes have driven L2 to be quite small. The 
second term in the right-hand side of the above equation therefore dominates, making 
P2/T a decreasing function of time. So, at the start of the low temperature epoch, 
P2/T slowly decreases, converting u^s and i/^'s to i/^'s as it does so. As for the Group 
1 modes above, it is easy to see that if the two Group 2 modes had slightly different 
resonance momenta, P2 and P|, then the dynamics forces ~^ ^2- The effect of the 
three-flavour z/g — t'^ — u'^ subsystem is to convert |z/^) to -^iWe) + I'^a*))- So, as P2 
moves (backward) through the neutrino momentum distribution. 
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(106) 



The effect of this conversion is to generate significant L,y^^ and L,^^ asymmetries which 
are negative in sign (given that we have taken L,y^ > for definiteness) . As the 
evolution unfolds, at some temperature P2/T changes direction and begins to increase 
again. This is due to the gradual increase in L2 which eventually makes the second 
term on the right-hand side of Eq. ( |105D smaller in magnitude than the first term. 
The full analysis, incorporating all of the modes simultaneously, requires a numerical 
treatment. We compute the minimum P2/T to be ~ 6, which is still in the tail 
of the distribution. The upshot of this somewhat complicated evolution is that the 
P2 resonance momentum does not sweep through the entire momentum distribution, 
but rather, it sweeps through a significant part of the high momentum tail. In the 
temperature regime where P2/T makes the return journey from its minimum value 
back to high values, the adiabatic MSW transitions have little effect because they 
simply swap the almost equal number densities of u'^ and i^e./x that were created by 
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adiabatic transitions before the turnaround.0 This is of course only true provided 
that the momentum distribution of v'^ states does not get significantly modified by the 
mirror weak interactions, an issue we will discuss in more detail later. 



The Vr ^ Ve O'lT'd Vr ^ Group 3 modes. These modes, being ordinary-ordinary, are 
slightly different in character to the ordinary-mirror modes. Their effect is to reprocess 
some of the L^^ into L^,^ and Ly^ . In the early stages of lepton number creation, all of 
the ordinary-ordinary modes are unimportant, because they simply swap flavours with 
almost identical number density distributions. However, eventually the L^^ asymmetry 
created by the Group 1 modes is large enough to distort the Vt momentum distribution 
so that the attendant reduction in the Vj. number density relative to that for Fg./i allows 
V-r ^ Ue^fj. oscillations to induce nontrivial dynamics: the depletion of and 17^ states 
at the resonance momentum P3. This effect becomes significant when L^,^ becomes 
quite large, which occurs roughly when Pi/T ~ 2. Now, from Eg. (pQ]) it is evident 
that P3 ~ |Pi using the fact that L^^ is the largest lepton number in the system. 
This has the important consequence that the overall effect of the Group 3 oscillations 
is not very large because by the time Li,^ is large, P3 is already well into the tail of 
the momentum distribution. This is fortunate, because these oscillations are more 
complicated to describe. Unlike the Group 1 and 2 modes, it is easy to see that 
if the resonance momenta, P3 and P3^, of the two modes are slightly different, then 
they do not subsequently evolve to coincide. (The reason is that if, say, the 17^ 



resonance momentum preceeds the 



resonance momentum then the u^- 



oscillations act to reduce L^'^^^ thereby increasing the rate at which this resonance 
momentum moves relative to the Vt ^ resonance momentum.) Note, however, 
that because of the more influential Group 2 oscillations, it follows that P3 and P| are 
at least approximately equal. What precisely happens will then depend on the width 
of these resonances and whether they overlap or not. This means that the effects will 
be dependent on the values of the relevant oscillation parameters. In our numerical 
work we will assume that the resonances overlap, so that IVj.) ^ TJ^I'^'^^ ~'~ I^m))' 
this case. 



NvAPs) 



^ + ^ + iV..(P3) 

Nv^iPs) , N^^XPs) 



+ 



2 2 

iV^^(P3)+iV. (P3) 
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We stress that, were the above assumption prove to be invalid, our numerical results 
would not be greatly affected because the Group 3 modes have a relatively weak effect 



Actually, in this region, the two resonance momenta P2 longer dynamically driven 
to coincide, making the oscillations somewhat more complicated. However, if the ly'^ tail is fully 
populated from the previous evolution of the system then this complication matters very little. 
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for the reasons discussed above. Finally we note that the L^^ and L^^ asymmetries 
created by the Group 3 modes have the opposite sign to the L^^ and L^^ asymmetries 
generated by the Group 2 modes. 

4. The V'^ <-i> 77g and V'^ V'^ Group 4 modes. These mirror-mirror modes can be ne- 
glected because P4 is always greater than Pi, and thus the u'^^^ states are approximately 
empty (as is V'^) when the P4 resonance momentum moves through. 

Having understood to some extent the effect of each group of oscillations, it is now time 
to solve the complete system of coupled equations for the various lepton numbers. These 
are obtained by a straightforward generalisation of the two-flavour case given in Eg. (PD|) . 
They are: 
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where 



Xo 



T 
T 
T 



N^APi)-\[N-,.SPi) + Np^iPi). 



X3 = - X^,(P3) - ^[NvM) + N^,{P,)] 



Expanding out Eq.( |lQ8| ) we find 

dL,, 



y2 



dT 

dL^^ 

~dT 
dLyi^ 
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a + p- 
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dLyi 



dT ' dT 

dLy.^ ^ dLyi 
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;iii) 



and /i = 1 for d{Pi/T)/dt > and /i = -1 for d{Pi/T)/dt < (i = 1, 2, 3). Solving Eq. ([lTOD 
we find, 





(^2/3 + Cv)iyiy3 - 70) + (p2/3 + 00 ("2/3 + iv) 




(2/22/3 -C^)(2/i2/3- 70) - 
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;ii2) 



We compute tlie number densities incorporating Eqs.( |103] - [L07| ). Tlie repopulation and tlier- 
malisation of tlie neutrino momentum distributions is taken into account using the same 
expression as in the high temperature epoch: 



dtN'^iip, T,0) 



r.(p) 



N'''i{p,T,0) Ar«i(p,r,0) 



(113) 



where T^ip) is the total collision rate. The chemical potentials are computed from the lepton 
numbers as per Sec. III. 

Observe that the oscillations will necessarily generate a significant number of v'^,v'^ 
and v'^ mirror neutrino species. We will make the simplifying assumption that there is 
negligible thermalisation of these mirror neutrinos. By this we mean that the mirror weak 
interactions of these mirror states are weak enough to not appreciably modify the mirror 
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neutrino momentum distributions. We will discuss later the circumstances required for this 
is to be a vahd approximation, and the expected effects when it is not vahd. 

In solving Eq.( |112| ) initial conditions for L^^ and L^i^ for each a must be specified at 
a temperature, Tiow, that serves as the initial point for the low temperature epoch. For 
definiteness we take Tiow = where Tc is the critical point during the high temperature 
epoch at which the explosive growth of L^^ begins. Our results are quite insensitive to the 
precise value taken for Tiow so long as it is high enough for the lepton numbers to be still 
much less than 1. (This issue is discussed more fully in Ref. p6[|.) The lepton numbers are 
related to the resonance momenta by using Eqs.(PBD and 



_ ^ / 2 1 9 

_ u; / 2 1 15^ 

''^ ~ 24 VP^ ^ ^ ^ 

u (U 17 15\ 

where u = ^^^^t ■ (The Group 4 oscillation modes have been neglected for reasons discussed 
earlier.) Thus, specifying the values Pi/T at T = Tiow completely fixes the values of the 
lepton numbers at that temperature. From our numerical work, we find that the T ~ Tiow 
values of the resonance momenta Pi and P2 are given approximately by 

Pi/T ~ 0.3, P2/T ~ 15. (115) 

These values are approximately independent of the vacuum oscillation parameters so long 
as the various mixing angles obey sin^ 26 ~ lO"^'' and provided the (5m^'s lie in the range of 
interest. Also note that our subsequent numerical work is not very sensitive to the precise 
initial values of Pi/T and P2/T provided that Pi/T is small (less than about 0.6) and P2/T 
is large (greater than about 10). We also need to specify the initial values of the signs fi. 
We take /i = /s = 1 and /2 = — 1 at T = Tiow- Subsequently fi are evaluated from the 
previous time step. 

In the region during and just after the exponential growth, the initial production of L^^ 
and due to the oscillations z7^ ^ Ve and 17^ ^ z7^ is suppressed because the number 
densities of all the ordinary neutrino flavours are almost equal. At T = Tiow we find that the 
creation of lepton number due to these oscillations is approximately negligible. This means 
that the main contribution to L^^, L^^ and L^/^ at T ~ Tiow is from u'^ ^ z/e,^ oscillations, 
and thus L,y'^ ~ —2Ly^. It follows that the initial value for P3/T can be approximately 
related to the initial values of Pi/T and P2/T by 

P3 33Pi 33P2 ^ ' 

We have solved this system of equations for the illustrative example of Sm^^^^^ = 50 eV^. 
In Fig. 3 we show the evolution of the four resonance momenta Pi/T. The evolution of L^^ 
and Lyi^ for the same ^^«farge parameter choice is plotted in Fig.4 for the Ly^ > case. 

Let us now turn to the implications of the oscillations for BBN. The change in Yp due 
to the neutrino oscillations can be separated into two contributions. 
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5Yp = 5iYp + 52Yp, (117) 

where 6iYp is the change due to the effect of the modified electron neutrino momentum 
distributions on the reaction rates, and 52Yp is due to the change in the energy density (or 
equivalently the change in the expansion rate of the universe). The former effect can be 
determined by numerically integrating the rate equations for the processes given in Eq.(p5D 
using the modified electron neutrino momentum distributions N^^ and N^:^ as discussed in 
Appendix A. The latter contribution can be computed from the momentum distributions of 
the ordinary and mirror neutrinos through 



5, 

where 



Yp ^ 0.012 r[KM + N^Av) + K'Sp) + N^'M]pdP ' ^ ' (^S) 



po ^ / N^'^ip, T, 0)pdp = —T\ (119) 
JO 240 

is the energy density of a Weyl fermion at equilibrium with zero chemical potential. [Recall 
that Eq.(|86D can be used to express 5Yp, 5iYp and 52Yp in terms of effective neutrino number, 
(^A^z/.off) ^iNy^eS and (52A''i.,efr, respectively.] To calculate (52^p, "we numerically determine the 
momentum distributions at T = 0.5 MeV. Because of the approximate kinetic decoupling 
of neutrinos for temperatures below about 3 — 4 MeV, large contribution^ to 52^p, should 
they exist, must have been generated earlier. A temperature of 0.5 MeV is therefore a safe 
place to evaluate the final 52Yp. 

Recall that there is an ambiguity concerning the sign of the Lj^^ lepton asymmetry. We 
have considered the Ly^ > case above for definiteness, but L^^ < is equally likely a 
priori. (See Refs. ||2^ for further discussion of this.) For the negative Lj^^ case, the roles 



of particles and anti-particles are reversed for the modes quoted in Eq.(|96|) and subsequent 
equations. One consequence of this is that the signs of all the other asymmetries are also 
reversed. The quantity 5iYp will obviously be significantly affected by this ambiguity in 
sign, while ^2^p will not be affected at all. This means that we have two possible values for 
the overall change in the effective number of neutrino flavours during BBN. The results of 
the numerical work is presented in Figs. 5 and 6. Figure 5 treats the L^j^ < case (which it 
turns out means that L^^ > 0), while Fig. 6 displays the L^^ > case (which implies that 
Lu. < 0). 

Observe that 6iYp is not very large. The main contribution to it is from the modification 
of the high momentum tail of the z/g distribution due to the Group 2 u^. ^ oscillations. 
This is partially offset by the modification of the distribution due to Ur ^ Ve oscillations. 
It is also evident that 52Yp is close to zero for 5ml^-^^^ ~ 300 eV^. This is simply because 
the generation of mirror states, which is dominated by the v'^ and Vr ^ ^'^ modes, 

occurs below the kinetic decoupling temperature for i/^'s. This means that the z/^ states 
which have converted into mirror states are not repopulated. For larger values of ^^«farge' 



^^By 'large contributions' we mean 52-^i/,efr ~ 0.10. 
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the Vr states begin to get repopulated, and the energy density increases accordingly. We 
should also emphasise that our calculations contain approximations. The most important 
are that the repopulation is handled approximately via Eg. ( |113| ) and we have neglected 
mirror thermalisation. Thus, our results have a theoretical uncertainty, which we estimate 
to be of order dN^^eS ~ 0.3 (see following discussion). 

We now discuss in more detail the effects of mirror neutrino thermalisation. Recall 
that in the foregoing computations, we have included mirror neutrino thermalisation via 
the quantity 7p = (T'/T)'* during the High Temperature Epoch, but neglected it during 
the Low Temperature Epoch. (Remember that mirror weak interaction rates increase with 
temperature.) We now discuss when this approximation is valid and the expected effects 
when it is not. 

Given our division of the evolution of the system into High and Low Temperature Epochs, 
it is convenient to also classify mirror weak interactions into two categories. The first cate- 
gory consists of the interactions of the mirror neutrinos generated during the Low Temper- 
ature Epoch with the background mirror neutrinos, electrons and positrons left over from 
the preceding High Temperature Epoch (henceforth designated MHTB for "mirror high-T 
background"). The second category is the elastic collisions of the 17^, v'^ and v'^ neutrinos 
generated during the Low Temperature Epoch with themselves. We now estimate each of 
these thermalisation rates. 

The interaction rate of a mirror neutrino v'^ of momentum p with the MHTB is approx- 
imately given by, 

12 T^S / P 



Now, from the discussion above [see Eq. (|96D ], the additional mirror neutrino states created 
during the Low Temperature Epoch consist of the flavours v'^, v'^ and z/^. Their interactions 
with the MHTB can be approximately neglected if 

^^^1 ^ l.y.Gl'^i^^T'^l. (121) 
H ^5.5 V3.15T/ ^ ' 

As summarised in Fig.3, the resonance momentum P2IT for the modes producing z/:J.'s is 
always higher than the resonance momentum Pi/T for the modes producing v'^ and v'^. The 
mirror thermalisation effects will therefore be most important for the v\ states. It is also 
clear that the relatively high momentum v'^ states are produced at a higher temperature 
than z7g and V'^ states of a corresponding momentum. The temperature range of interest 
for the v'^ lies between Tiow — T^l and the temperature Tmin at which P2IT reaches its 
minimum value of ~ 6. 

We can easily numerically compute Tmin to obtain, 

iz!!l^L~0.70 . (122) 

MeV V y ^ ' 

We now estimate the effects of thermalisation by considering the interaction rate for a v'^ 
with a typical momentum of p/T ~ 8 at a temperature around or slightly higher than Tmin- 
Although the constraint Eq.( |121| ) on 7p is stronger for higher values of T, the number of 
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z/^'s produced is lower, so the effects of their thermahsation will be correspondingly weaker. 
The choices made for p/T and T as input for Eg . (11211) represent a reasonable "compromise" 
driven by these considerations. So, we estimate that the interactions of v'^ with the MHTB 
can be neglected provided that 

< 1 / eVM^ , , 

We now estimate the thermahsation rate, due to elastic collisions with themselves and 
with the z7g^^ states produced during the Low Temperature Epoch, of the mirror u'^ states 
produced during the Low Temperature Epoch. We will collectively call the mirror neu- 
trino/antineutrino states produced during the Low Temperature Epoch as the "mirror low-T 
background" (MLTB). Let us denote the collision rate for a v'^ of momentum p with the v'^ ^ 
{v'^) component of the MLTB by Ti{p) [r2(p)]. The relevant collision rates can be obtained 
from Ref. ||4^ to yield 



Ti{p) ^ O.ISG^T^ 

r,(p).0.77G^T=(^;lJ^^J, (124) 

where po is defined in Eg. ( |119D . So, the interactions of the v'^ with the 17^^^ (i/^) MLTB can 
be approximately neglected provided that 

^^i^l ^ 0.13G^T3^f^^^V^Ui, 
H ^ 5.5 \ po J \3.15Tj 



'Pv' + Pv\ 




H ^ 5.5 V po / v3.i5ry 



(125) 



< 



Our numerical work shows that {p-pi^ + p^'^)/ <^ 1 until guite low temperatures T/MeV 
(^mfj^rge/sV^)^/^. In fact the second condition in Eg. ( |125| ), from u'^v'^ elastic collisions, 
is the more stringent reguirement. For this case we can estimate the collision rate by 
considering a v'^ of typical momentum p/T ~ 8, which is produced at a temperature 
T/MeV ~ 1.3(5m2^,g,/eV2)V4. The ratio of energy densities reguired is estimated from 
adiabatic conversion as P2/T evolves from its initial value to about 8: 

p< - r y'^^y ~ o.o4. (126) 



po 27r^po ^8 I + ey 



Using these numbers we estimate from Eg.( |125| ) that v'^v'^ elastic collisions can be approxi- 
mately neglected provided that 



0.05 -1 ^ 5mt^^^5QeY\ (127) 
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It turns out that this numerical bound is not very sensitive to what we choose for a "typical" 
P2/T, so it is fairly robust. Thus, in summary, we conclude that the thermalisation of the 
mirror neutrinos can be approximately neglected provided that 



3 



\J<^J ^"""^ 5<a.ge-50eV^ (128) 

Let us now discuss what happens when there is significant thermalisation of the mirror 
neutrinos. Let us first consider the case of the v'^ states. If they are thermalised, the v'^ 
distribution will be close to an equilibrium distribution given by 

1 p2 

Nl} = ^. (129) 

The elastic u'^ collisions with the background v'^ conserve both the number density n^,^ and 
energy density p^,/. , so these quantities can be used to determine the two parameters /z^,/. and 
7;;. Now, 



/ iV.;prfp^^e"<. (130) 



Thus, 



Hence for any T, we can compute T^,/. and [i^i^ by computing p^i^ and nyi^ from N^i_^. To 
estimate the effects of the repopulation we use the equation. 



dt N'^'iip, T, 0) 



repop 



N'''i{p,T,0) A^^q(p,T,0) 



(132) 



where T2{p) is the elastic u'^u'^ collision rate quoted earlier. 

The most important effect of the thermalisation of the u'j. occurs when P2/T evolves 
from P2/^|min to infinity. For the example of Figs. 3 and 4, namely 5^fargc ~ eV'^, 
we have computed A^^? at the temperature T^in- This is shown in Fig. 7 as a function of 
p/T va. order to compare this distribution with the distribution of u'^ states which would 
exist in the absence of v'^ thermalisation. In this latter case, the MSW transitions that 
occurred during the previous evolution of P2/T from about 15 down to about 6, populated 
the v'^ states from the tail of the z/g.^t distributions (which have approximately negligible 
chemical potentials in this region). Furthermore, ordinary weak interactions repopulated 
the depleted z/g,/* tails. This means, in the absence of mirror thermalisation, the journey 
back from P2/7'|min to infinity is dynamically inert as the oscillating species always have 
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approximately equal number densities in the resonance region. However, if the i/^'s are 
thermalised this is not the case. 

Computing the evolution of the system when P2/T evolves from P2/T\min back to high 
values appears to be problematic. The problem is that in this region the resonance momenta 
P2 and P2 for the modes v'^ ^ Ve and v'^ ^ v^, respectively, are not dynamically driven to 
coincide. Thus, in this case, we might expect different results depending on which resonance 
momentum goes first. Since the previous evolution of the system was such that the two 
resonance momenta coincided, it is not clear which resonance momentum will in fact go 
first. For instance, the result may well depend on a statistical fluctuation, and therefore 
may be different in different regions of the universe. The physical implication of this would 
be a spatially dependent Yp distribution. 

We have made some numerical estimates using the prescription given in Eq.( |106| ). Our 
numerical results indicate that the overall affect of v'^ thermalisation is not unacceptably 
large, typically about ^A^'j.efr ~ 0.3 [for the entire range of interest in Smf^^g^. For L^^ > 
the effect is positive, that is SN^^cS ~ +0.3, while for L^^ < the effect is negative, that 
is SNjy^cS ~ —0.3. This essentially results in a theoretical error of this magnitude for the 
parameter space region which violates the bounds in Eq. ( |1 28| ) . 



C. Low temperature neutrino asymmetry evolution in the EPM: Case 2 



We now consider the case where the z/^ and u'^ masses are not negligible. This is of 
considerable interest since mu^^ ~ 1 eV is expected if the LSND anomaly is due to 
neutrino oscillations. We consider the mass hierarchy, 

m„^+ - m^^. > "^!.^+ - > my^+,m^,_- (133) 

In Case 2, there are two 5m^ scales. The modes listed in Eq.(P3D have the large 5m^ = 
5mf^,g„ while 

J^t^ ^ J^t^ ^ J^e, Jy'f, ^ Z^e, Jy'^ ^ V'e, (134) 

plus the associated antiparticle modes have the smaller 5rn? = Sm'^^^ii. The specification of 
the parameter space of interest is completed by taking the remaining modes, ^ u'^ for 
a = e,fi, T, to have a negligible (much smaller than 5^«smaii)- 

The four modes in Eq.( |134| ) typically have distinct resonance momenta which we denote 
as follows: 
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= Pl, 




<— > 
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= P2, 






Ve, 
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= P3, 






— / 


Pres 


= Pa- 



(135) 

The resonance momenta for each of these modes can be obtained from Eq.(p7D, 

Pi _ ^"^small (136) 
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where z = 1, . . . , 4 and 



(137) 



Observe that we use the lower case pi notation for the Sml^^^^ modes of Eq.( |135| ) and the 
uppercase Pi notation for the 5mf^^^^ modes of Eg . (p6|) . 

Since 5^«smaii ^ ^"^large; foUows that pi^2 ^ Pi- (Note however that ps and p4 start 
out being infinitely large because of the L^^ = L^^ and Lj,/ = L^/^ conditions. As we will 
explain shortly, these two modes have little effect). For our numerical work, we will consider 
the parameter space region where the hierarchy between ^rn^j^^^^ and Sniig^j.^^ is great enough 

so that Pi,2/T ~ 0.5 when Pi/T ~ 10. Numerically, this corresponds to Smf^^^^ ~ SO^m^jj^g^y. 
Let us denote by T = the temperature at which the Pi/T are all greater than 10. From 
Fig. 3 and Eq.(|97|) it is easy to see that is given by 



The evolution of the neutrino ensemble during the Low Temperature Epoch for Case 2 
therefore breaks up into two temperature regions: T ~ and T ~ T^. When T ~ T^, 
the evolution of the system is dominated by the ^^«iarge modes of Eg . (|96|) , because the 
<^^smaii modes are negligible due to their very small resonance momenta. In this temperature 
region the evolution of the lepton numbers and number densities can be evaluated using the 
equations of the previous section. When T ~ T^, the 5w^farge modes are no longer effective, 
because they all have Pi/T ~ 10. The ^w^smaii modes begin to become important. For 
~ 1 eV^, it follows from Eq. (|T36|) that Pi,2/T ~ 1 for T ~ 1 MeV. This means that 
for 5/«smaii m the LSND range, the oscillation modes of Eq.( [134| ) become important while 
the BBN reactions n ^ p are still rapid. So, these oscillations can potentially influence 
BBN and therefore should not be ignored. 

It turns out that it is not possible to use the adiabatic approximation [as encoded in 
Eq.(PDD] to work out the effects of the oscillation subsystem given in Eq. ( |134| ). This is 
because of the structure of Eq. (|136|) . For example, z/^ z/g oscillations create significant 
L^^ asymmetry, as P2/T sweeps through the momentum distribution. The L^^ asymmetry 
becomes so large that L2 0. This makes the rate of change of P2/T very large and the 
system is no longer adiabatic. Because of this complication, we will analyse the effects of 
the modes in Eq.( |134| ) using the Quantum Kinetic Equations. 

We start integrating the QKEs for the subsystem of Eq.( [134| ) at T = Tj. with the values 
of the number distributions A^^^ and N^'^, and the lepton numbers L^^ and L^'^ (a = e, fi, r) 
obtained from the evolution equations of the previous Subsection. We will not explicitly write 
down the QKEs for this subsystem here, because their form is obvious once the contents of 
Sec. IV above have been understood. Nevertheless for completeness we will include them in 
Appendix B. The reader, however, should not the following points: 




(138) 
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1. We utilise the approximation that the x and y components of the polarisation vectors 
for the modes of Eq. (|134|) vanish at T = T^. While this is not expected to actually 
be the case, the subsequent evolution is not sensitive to the particular choices made 
for the initial values. This is because the resonance momenta for the subsystems are 
either very small or very large, so the very first stage of the evolution after T = 

is fairly unimportant. Furthermore, correct values for the x and y components are 
quickly generated by the QKEs soon after T = T^. 

2. Repopulation and thermalisation of the mirror neutrino ensembles have been neglected, 
consistent with our treatment of the evolution between the end of the High Tempera- 
ture Epoch and T = T^. Actually, this is an excellent approximation in this regime, 
since for typical interesting parameter choices the mirror sector temperature is quite 
low. 

3. The ordinary-ordinary 17^ and mirror- mirror u'^ u'^ modes can, to a good 
approximation, actually be omitted. Recall that the resonance momenta ps and 
are initially very large, much larger than pi and p2 respectively. Subsequent evolution 
maintains this hierarchy in the resonance momenta. The z7^ ^ 77^ mode, with reso- 
nance momentum pi, is strongly reprocessing lepton number as the resonance moves 
through the body of the distribution. The coupled ^ Vg mode, on the other 
hand, sees its resonance momentum p^ remain in the tail of the distribution, where 
it is ineffective because of the essentially identical number densities of and Vf. in 
the tail. Remember that because weak interaction rates after typical values of are 
getting quite weak, there is little thermalisation of the reprocessed e and n asymme- 
tries created at low momenta pi. In other words, the p^ resonance barely "knows" the 
asymmetry is there. 

In Figs. 8, 9 we plot the evolution of the resonance momenta, Pi/T and Pi/T for an illus- 
trative example. We choose 5w^iargc = 50 eV^, 5w^smaii = 1 ^V^ and all of the vacuum mixing 
angles to be 10^^. We emphasise that our results should be approximately independent of 
the vacuum mixing angles so long as 10^^° ~ sin^ 2^ ^ 1. In Fig. 10 we plot the evolution 
of all of the asymmetries for the same example. 

The effect of the oscillations on BBN is given in Figs. 11- 16. In Fig. 11 we have plotted 
Ny^es versus Sml^^]^, with 5mfj^j.gg = 50 eV^ and L^^ < 0. Figures 12 and 13 are similar except 
'^"^large = 200 and 800 eV^, respectively. Figures 14-16 are the same as Figs. 11-13 except 
that the opposite sign asymmetries have been considered. It is very important to note that 
the effect of a nonzero Sm^^^n is considerable. In particular, as Figs. 11-13 illustrate, N^^^ff 
depends sensitively on Sm^^^n, with negative corrections to N^^^ff equal to about one effective 
neutrino flavour are achieved for Sml^g^n values in the few eV^ range. In Figs. 14-16 observe 
that has a large effect even at very low values compared to the corresponding effect in 
Figs. 11-13. This asymmetry is due to the neutron/proton mass difference. 

Recall that Figs. 11-16 have not included the effects of mirror neutrino thermalisation. 
As already discussed, these effects should be significant for the modes with = Smf^^^^ 

if ^^Ti^arge ~ ^0 cV^. The thermalisation of mirror neutrinos should have negligible effect for 
the modes with (5m^ = Sm'^^g^i because they are only important when the temperature is 
typically less than about 1 MeV. As discussed earlier, our rough estimate of the effect of 
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the mirror thermalisation is about SN^^cS ~ 0.3. For L,y^ < {L,y^ > 0), the effect of mirror 
thermahsation should be to decrease (increase) SNy^eS by of order —0.3 (+0.3). 

As mentioned above, the numerical results were obtained using the parameter choice 
sin^ 29^e' = 10~®. Actually we expect the results to be quite insensitive to sin^ 29^^' so long 
as sin^ ^ 1- The reason is that the amount of L^^ that gets created is already close to 
the maximal amount possible. That is, after its rapid creation (which is at T ~ 0.6 MeV 
in the example in Figure 8), the quantity L^^f^^ ~ 0. Increasing sin^ 2^^e' cannot increase 
the amount of L^,^ much since it is already close to the maximum possible. Also, it cannot 
be created much earlier. Thus, the results shown in Figs. 11- 16 should be approximately 
independent of sin^ -0 

Finally, we should remark that the results of this section indicate that the bounds ob- 
tained in Sec. IV can be evaded somewhat. The reason is that even if 6N,^^ea — 1-5 from 
the high temperature population of mirror states from the I'r — t'^ oscillations, this can 
be compensated by a SNi,^^^ ~ —1.0 from the low temperature generation of a large L^^ 
asymmetry. 



VI. IMPLICATIONS FOR HOT DARK MATTER 

In the scenario considered in this paper, where v'^ oscillations solve the atmospheric 

neutrino anomaly, a BBN bound of Ny^^s ~ 3.6 implies ~ 1 eV for l^m^^-jj^^J ~ 10~^'^ eV^ 
(see Fig. 2). Neutrino masses in the eV range have long been considered cosmologically 
interesting, because they would make a significant contribution to the energy density of the 
universe. In the standard Big Bang model, the contribution of massive standard neutrinos 
to the energy density is given by the well known formula, 

/i292 eV' ^ ^ 

where h is the usual cosmological parameter parameterising the uncertainty in the Hubble 
constant. Thus, neutrinos in the eV mass range are a well known and well motivated 
candidate for hot dark matter. 



Before the advent of information at high red shift values [50|, large scale structure for- 



mation studies strongly favoured a hot plus cold dark matter mixture with VLy ~ 0.20 — 0.25 



|51|. While recent work incorporating the new high redshift large scale structure data has 
reduced the need for a hot dark matter component, it remains an interesting possibility. 
Given that Vr masses greater than a few eV or so are well motivated from the combined 
requirements of the atmospheric neutrino anomaly and BBN, we see that the existence of 
neutrino hot dark matter is a generic prediction of the EPM. 

There is an interesting complication in the hot dark matter story for the EPM (and 
models with sterile neutrinos) which we now discuss. We will take by way of concrete 



We have numerically checked this by looking at the case sin^ 20^e' = 10 ^ and found almost 
identical results. We have also checked smaller sin^20^e'- For sin^ 20^e' ~ 10~^ the oscillations 
begin to become so non-adiabatic that the oscillations start to become less effective. 
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example that only the v^- (and z/^) has an eV scale mass. Again for the sake of the example, 
we will consider the neutrino mass range required by what was the favoured hot plus cold 



dark matter scenario fpl 



3 eV ~ m^^ ~ 7 eV, (140) 

even though the present situation is less clear. The point we want to make is that whatever 
a "favoured neutrino hot dark matter mass range" might be at any given time, the situation 
is modified somewhat in the case of the EPM model. The reason is that the ^ u'^ and 
Ut- z/g oscillations generate such a large Li^^ that the total number of tau neutrinos is 
actually significantly reduced. For the parameter region 

10 ~ 5ml^^JeV^ ~ 300, (141) 

the final value of L^^ is about 0.27. The large final lepton number occurs because about 
70% of the anti-neutrinos have been depleted (for the L^^ > case) which means that 
the total number of tau neutrinos plus tau antineutrinos is roughly 0.65 of the standard 
expectation. (Note that the total number of neutrinos has not changed much: the missing 
heavy tau antineutrinos have just been converted into light mirror states.) Also note that 
a small number of i'!^ are also generated by the oscillations, and it turns out that the total 
number of Ur and u'^ (plus antiparticle) states is about 0.70 of the standard expectation. 
The effect of this is to change the "favoured hot dark matter mass range" from what the 
expectation would be in the absence of mirror (or sterile) neutrinos. We can guess that in 
the context of the EPM the "favoured" tau neutrino mass is actually about 50% larger than 
the naive expectation. (A full large scale structure computation would need to be performed 
to fully explore the consequences of a depleted z/^ distribution.) Thus, in the EPM model, 
the hypothetical favoured mass range of Eq.( |140| ) becomes instead 



5 eV ~ m^^^ ~ 10 eV. (142) 
This means that Smi^^g^ is expected to be in the range 

25 eV^ ~ Sml^^^ ~ 100 eV^ (143) 

provided of course that the scenario of Eq.( |140| ) is correct. The hypothetical hot dark 
matter region of Eq. (|143|) is the shaded band on Fig.2. From this Figure, we see that there 
is considerable overlap between the BBN allowed region and the hot dark matter region. 

Finally, note that structure formation outcomes in hot plus cold dark matter models are 
generically sensitive to the number of eV neutrino flavours, not just to fl^- These studies 
typically assume that the number of eV neutrino flavours (usually taken to be degenerate in 
mass) at the epoch of matter-radiation equality is an integer. It is important to understand 
that this is only true provided that mirror or sterile neutrinos do not exist. Indeed, as we 
have just explained above, we expect N^'^^'^y ~ 0.70 in the EPM in the parameter space 
region of Eqs.(|l],|2D and ([Till) . 
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VII. IMPLICATIONS FOR THE COSMIC MICROWAVE BACKGROUND 



During the next decade or so, high precision measurements of the anisotropy of the 
cosmic microwave background (CMB) will be performed by several experiments (such as the 
PLANCK and MAP missions). These satellites should be able to measure detailed spectral 
properties of the electromagnetic radiation in the universe at the epoch of photon-matter 
decoupling ^2 . In this context it is important to note that mirror and sterile neutrinos 



can leave their 'imprint' on the cosmic microwave background [Q. This information will 
complement knowledge obtained from BBN because, (i) BEN and photon decoupling take 
place at different epochs, and (ii) BBN is sensitive to both the expansion rate and the direct 
effect of L^^ on nuclear reaction rates whereas the CMB is insensitive to the direct effects 
of the asymmetry. Because of point (ii) we have to distinguish between expansion rate and 
effective neutrino flavour counting. So, it is useful to introduce the quantities iV^'^^* and 
jyheavy -y^^j^jdi effectively count the number of light neutrino and heavy neutrino flavours, 
respectively, at the epoch of photon decoupling. These quantities, which quantify expansion 
rates, are to be used in conjunction with Nj^^^s which contains both expansion rate and L^^ 
information. It is important to appreciate that the number of relativistic neutrino flavours 
may be different at the time of photon decoupling compared to BBN. So, in this context, 
'light' means much less than about an eV, making these neutrinos relativistic at the epoch 
of photon decoupling, and 'heavy' means more than about an eV, making those neutrinos 
approximately non-relativistic. Of course in the minimal standard model of particle physics 
with its three massless neutrinos, A^j^^efr = N^^^^ = 3 and N^'^^^^ = 0. However, in models 
with sterile or mirror neutrinos, N,y^eS 7^ A^^'^^* and N^'^^^^ 7^ in general. It is also important 
to appreciate that in the EPM (or in models with sterile neutrinos), none of these quantities 
is in general an integer. 

The CMB implications of the EPM depend on the neutrino parameter region. If we take 
by way of example the mass hierarchy of Eqs. (|9l[) and (^2]), with rriy^^ ~ 1 eV (as suggested 
by Fig. 2) then 

heavy ^ + n,, + autiparticles .ght _ P^e + Pu, + P.^ + Pw^ + antiparticles 

2no ' ^ 2po ' ^ ' 

where rii {pi) is the mass (energy) density of species i with uq (po) being the mass (energy) 
density of a Weyl fermion distribution with zero chemical potential. Taking ^m^^rge in the 
range Eq. ( |141| ) we find that 

^heavy _ q ^q^ flight _ 2.3. (145) 

This should be distinguishable from the minimal standard model expectation. 

We conclude by emphasising that in general the precise measurements of the CMB may 
well prove to be quite useful in distinguishing between various competing explanations of 
the neutrino anomalies, since each model should leave quite a distinctive imprint on the 
CMB. 
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VIII. CONCLUSION 



The Exact Parity Model is theoretically well motivated by the neurotic desire of some to 
have the full Lorentz Group as an exact symmetry of nature. It is very interesting that this 
model can, essentially as a byproduct, provide an elegant explanation of the atmospheric 
and solar neutrino problems in a way that is fully compatible with the LSND results. In this 
paper, we explored the novel cosmological phenomena implied by the existence of mirror 
neutrinos. 

We fo cussed on the parameter space region 

1^1^,+ ^ rn^,- ~ fn^^+ ^ fn^^,- ~ rn^r+ - ^y-r- (146) 

with all intergenerational vacuum mixing angles obeying 

10^1° ~ sin2 20 < 1. (147) 

The mass splittings amongst the e-like and /i-like states were chosen to solve the solar 
and atmospheric neutrino problems, respectively. The evolution of the neutrino and mirror 
neutrino ensembles was then calculated for the cosmological epoch between T = and Big 
Bang Nucleosynthesis. Generic outcomes were obtained for significant regions of parameter 
space because (i) some of the final neutrino asymmetries turned out to be independent of 
the oscillation parameters for a range of those parameters, and (ii) many of the modes were 
adiabatic and hence independent of vacuum mixing angles. 
The most important specific conclusions were: 

1. The — i> z/^ solution to the atmospheric neutrino problem is consistent with Big 
Bang Nucleosynthesis for the parameter space region illustrated in Fig. 2. The Vt 
mass implied by this region makes the Vr a hot dark matter particle. This calculation 
improves on that discussed in Ref. through the use of Quantum Kinetic Equations. 



The effect of EPM neutrino oscillations on the primordial Helium abundance has been 
computed. We find that a large change to the effective number of neutrino flavours 
during Big Bang Nucleosynthesis is produced for a range of parameters. In particular, 
a change equivalent to adding or removing about one neutrino flavour is obtained when 
the — ^fi mass splitting is in the LSND range. 
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APPENDIX A: DETAILS OF THE HELIUM ABUNDANCE COMPUTATION 



The modification of the v^. and z/g distributions due to the creation of Ly^ affects Big 
Bang Nucleosynthesis. This is primarily due to the modification of the n ^ p reaction 
rates. The result of this is a modification of the neutron/proton ratio. The most important 
observable effect of a small change to the neutron/proton ratio is a modification to the 
prediction for the Helium mass fraction Yp. This effect can be expressed as a change in 
the predicted N^^^s through the well known relation 6Yp ~ 0.012(5Ai'i,efr-^'' In computing the 
modification of Yp due to the modified neutrino distributions, N^^ and N-^^, we do not need 
to use a full nucleosynthesis code. The reason is that the effects of the modified neutrino 
distributions are only important for temperatures T ~ 0.4 MeV, well before nucleosynthesis 
actually occurs. A review of standard Helium synthesis which we found useful was Ref. . 
Our approach and notation follows this treatment quite closely. 

Recall that the primordial Helium mass fraction, Yp, is related to the ratio of neutrons 
to nucleons, by Yp = 2X„ just before nucleosynthesis. X„ is governed by the differential 
equation 

rl X 

^ = A(n-.p)X„-Ab^n)(l-X„), (Al) 



dt 



where 



A(n ^ p) = \{n + z/g — > p + e ) + \{n + e"^ — ^ p + z/g) + \{n p + e + z/g), 

X{p ^ n) = X{p + e~ ^ n + Ue) + X{p + Vf, ^ n + e^) + X{p ^ n + + Ue) ■ (A2) 

The rates for these processes are given by 

X{n + Ve^p + e~)= A [ dp^, 

J 1 + exp(— _£/e/T) 

l + exp(Ee/T)' 

l + exp(-Ee/T)' 

l + exp(^e/T)' 



A(n + e+^p + z7,) =A / ,^":L dpe, 



Xin-.p + e~+V,)=A I ^^^^^-^^,dpy. 



X{p + e- ^n + Ue)=A I -^PL^A^dp, 



Of course we are not saying that this equivalence is exact. It is not. The change in Yp due to the 
modification of the Uf> and Vf, distributions cannot be exactly represented as a change in A'^^^eff- This 
is because these two effects will have different impacts on the other primordial element abundances. 
However, because a small modification in the v^. and V(. distributions, or a small change in N^^^.^, 
primarily affects Yp, our use of the relation 5Yp ~ 0.0125A^jy^efr is reasonable. We prefer to express 
our results in terms of dNj^^^s rather than 5Yp just because 5Ni^^c,s is a more familiar unit. 
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J l + exD(—Ep T) 



+ exp{-EjT) 

X{p + e'+V,^n)=A [ , l'^' !t^)^, dp,. (A3) 

where Ve = Pe/ is the velocity of the electron (we use h = c = 1 throughout) and Nv 

is related to the neutrino distribution functions by Nv= 2tt'^N^. The constant A can be 
expressed in terms of the vector and axial vector coupling constants of the nucleon [J6], 



A = (A4) 



Also, Ee and Ei, are related by 





-E, 


= Q 


for 


n + z/g p + e~, 


E, 


-E, 


= Q 


for 


n + ^ p + T7e, 


E. 


+ E, 


= Q 


for 


n p + e~ +17^, 



(A5) 

where Q = rrin — rup ~ 1.293 MeV. The integrals of Eg. (|A5D are taken over all positive 



values of p,^ and Pe allowed by these relations. 

In order to compute Yp we need to know the time when nucleosynthesis occurs and 
neutron decay ceases. This is handled approximately by simply stopping the evolution of 
Xn at a point where agreement with the expected value of Yp ~ 0.24 occurs, which we find 
to be roughly when t ^ 300 seconds. This approximation does not affect the accuracy of our 
results at all since we are only interested in the difference between Yp using the modified 
and distributions and Yp using the standard distributions (i.e. Fermi-Dirac distributions 
with zero chemical potentials). Thus, to a excellent approximation, the modification of Yp 
due to the non-standard neutrino distributions has the form 

5Yp ^ 25Xn{t = 300s), (A6) 

where SXn{t = 300s) is the difference between X„(t = 300s) computed using the neutrino 
momentum distributions N^,^ and Njy^ and X„(t = 300s) using the standard momentum 
distributions. Of course the distributions N^^ and N-^^ typically depend on the time, so 
that the evolution of X^^ (t) must be computed concurrently with the evolution of N^^ and 
N-jy^. In solving the differential equation Eq.(AT), we employ the usual the initial condition 
xj= 0.5. 



We have checked our code against some previous calculations. For example, in Ref. p/7 
they consider the case of a time independent neutrino chemical potential (taken to arise 
from some unknown physics at high temperature). From Figure 2 of Ref. they find that 



6Yp ~ —0.020 for fi^/T = —fiu/T ~ 0.09 (for constant rj). Our code also gives exactly the 
same results under the same conditions. 

Finally note that at low temperatures T ~ mg, the e~^e~ annihilation process increases 
the temperature of the photons relative to the neutrinos. It also affects the time-temperature 
relation. In our numerical work, we take these effects into account using the equations given 
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in Ref. |^ (suitably modified to incorporate three light neutrino flavours instead of two). 
Of course this detail actually does not affect our results much, since most of the effects of 
neutrino asymmetries are only important for temperatures T ~ mg. Nevertheless, following 
E. Hillary (private communication) we include it because it is there. 



APPENDIX B: THE QUANTUM KINETIC EQUATIONS FOR THE MODES 



WITH 



JMImall in case 2 OF SECTION V 



This Appendix deals with the case defined by Eq.( |133| ), where there are a class of modes 



having Sm'fg^j.^^ [see Eq.(|93|)] and another class of modes having Sm^^g^^i [see Eq.( [134| )]. As 
discussed in Sec.V, the 5w^^arge ^"^smaii diodes approximately decouple from each other 
provided that Smf^j.^^ ~ bOSm^^^^^. The evolution of the ^rrii^j.^^ oscillations can be evalu- 
ated using the adiabatic formalism of Sec.VB. The 5?7ismaii modes can be neglected initially 
because their resonance momenta satisfy pi,2/T ^1. By the time T = [see Eq.( p.38| )], 
the 5"^smaii modes begin to be important, and provided that ^w^farge ~ ^^^"^smaii is satis- 
fied, the ^^rafarge modes can be neglected because Pi/T ~ 10. To compute the effects of the 
^'^smaii modes we must numerically integrate the quantum kinetic equations. Thus, we start 
the quantum kinetic equations at T = with the initial values of N^^, N,^/^, L^^ and Lyi^ 
obtained from the previous evolution involving the 5m\^^^^ modes. 

In this Appendix we do not follow exactly the notation of Sec. IV. We adopt an equivalent 
but simplifing change of variables which is very useful for complicated coupled-mode systems 
such as the one we are currently dealing with. 



For each of the four oscillation modes, we assign a density matrix 



In solving this system, it is convenient to use the variables P^, Py, N^^, Nj^^, N^/^ and N-p'^ 
rather than the variables Pq,P^, Py and P^. The Pq ^ are related to the A^'s as follows: 



Poip) ^ 
PoiP) 

Po'iP) 



n^M + Nk(p) 

A^<=<i(p,T,0) ' 

NM + N.'^{P) 
A^<=q(p,T,0) ' 

iV^.(p) + iV^.(p) 
A^<^q(p,T,0) 



P'ip) 
P'ip) 

pKp) 

P^P) 



iV^q(p,T,0) ' 

NM-N.'^P) 
A^°q(p,T,0) ' 

N'"i{p,T,0) ' 
A^^q(p,T,0) 



(El) 



The evolution of the number densities has a contribution from coherent effects and a con- 
tribution from the subsequent repopulation. Thus, 



dt 



dN,, 



dt 



+ 



dN„ 



dt 



repop 



dt 



dN^a 



dt 



+ 



dNjy 



dt 



(B2) 



repop 



The contribution from coherent effects can be broken up among the four modes as follows: 
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dt 

dt 

dN^ 
dt 

dN^, 
dt 



dN^ 
dt 



0, 



0, 

dN^ 
dt 



dNv, 
dt 

dN^: 
dt 



dN^ 
dt 

c 

_ dN^ 
dt 

dN^,, 
dt 

dN^> 
dt 



dN^ 
dt 



+ 



+ 



dN^ 
dt 

dNj,, 
dt 



dNj,,^ 
~dt 



(B3) 



with 



dN^,. 



dt 

dN^, 



dt 
dNv. 



dt 

dNr,' 



dt 



dN^' 



dt 

dN^ 
dt 

dN^^ 



dt 

dN^^ 



dt 



-^Jp,PlN^^{p,T,0)dp, 
~Jp2P;N^''{p,T,0)dp, 



y PsPlN^^{p,T,0)dp, 
^Jp,PlN^'^{p,T,0)dp, 



(B4) 



where Pi = —5m'^^^^YYsm26i/2p for i = 1,3,4 and = ^rngj^j^n sin 2^^j/2p for i = 2. The 
rate of change of the number densities due to repopulation is handled approximately via the 
equation, 



d N^Ap) 



dt N^'iip, T, 0) 



repop 



N^^^ip^T^O) N^^{p,T,0) 



(B5) 



where rc^(p) is the total collision rate, and N^{p, T, jj) is the equihbrium distribution which 
is a function of the chemical potentials which can be computed from the lepton numbers 
(see Sec. Ill for further discussion). 

The rates of change of the lepton numbers are given by 



dL„ dL^ 



dt 


dt 


dLu^ _ 


dL^^ 


dt 


dt 


dLi,i^ 


dLiyi 

e 


dt 


dt 


dL^i 


dL^' 



+ 



+ 



dt 



dt 



+ 



dK, 




dt 




dL^^ 




dt 








dt 




dLyi 




dt 





(B6) 
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where 



dt 
dLi,i 



dt 



dt 
dLyi 



dt 



dL^i 



dt 

dLu, 



dt 

dL^, 



dt 



dt 

dLyf 



= I PiPlN^''{p,T,Q)dp, 

^-^l(32P;N^''{p,T,0)dp, 



^ ' (3sPlN'''{p,T,0)dp, 



2n. 



I P,PlN'^\p,T,0)dp. 



(B7) 



Actually it turns out that the effect of the u'^ < 

nth 6171^ 



u' mode can be neglected because N-^' ipA ~ 
Nj;'^{p4). This is because the modes with 6171^ = Srn^^j.^^ create approximately equal numbers 
of u'^ and u'^ states. Also, the 77^ u'^ modes always have a lower resonance momentum 
I7g modes. This means that the change in N-p'^{p) due to 17^ 



u' modes 



than the 17^ <-> u'^ modes. This means that the change in N-p'^{p) due to 17^ v'^ 
does not occur until the v'^ <-> 1/^ resonance momentum has already passed by. By similar 
reasoning, the <-> Ve modes can also neglected to a good approximation. 
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Figure Captions 



Figure 1: \Ly^\/h (where h = T^/T^) versus temperature for ^ u'^^ oscillations with 
8w? = -50 eV^ and sin^ 29 = 10"^. 

Figure 2: Region of parameter space in the sin^2^T-^/, —Sm^^, plane where L,^^ is generated 
rapidly enough so that the z/^ ^ u'^ oscillations cannot significantly populate the z/^ states 

(for T ~ 0.4 MeV). This region, which in the figure is denoted by the "Allowed Region", 
includes all of the parameter space above the solid line(s). The top, middle and bottom 
solid lines correspond to the atmospheric Sm^ values of (5m^^,/eV^ = 10~^, 10~^'^ and 10~^ 
respectively. The dashed-dotted line is the nucleosynthesis bound Eg. (|i9|) , which takes 
SN,y^cS ~ 0.6 for definiteness, and the shaded region is the hot dark matter region indicted 
from some studies of structure formation (see Sec. VI). 

Figure 3: Evolution of the resonance momenta Pi/T, for the example with ^mfg^^gg = 50 eV^. 
The solid line, long-dashed line, short dashed line, and dashed-dotted line correspond to 
Pi/T, P2/T, Ps/T and P^/T respectively. 

Figure 4: Evolution of the lepton numbers for the same example as Fig. 3. 

Figure 5: SNi,^eS versus ^^«iarge Case 1 [see Eqs.( pl| , p^ )]. The dashed line is the contri- 
bution SiN,^^eS due to the effects of the L^^ asymmetry while the dashed-dotted line is the 
contribution S2Ni^^cS due to the change in the expansion rate. The solid line is the total 
contribution ^lA'^i.^efr + '^2A^i/,efr- This figure considers the case L^^ < 0. 

Figure 6: Same as Figure 5 except L^^ > is considered. 

Figure 7: A^^? (bottom solid line) at the temperature T = Tmin (see text) for the example 
of Figs. 3 and 4. The top solid line is the expected distribution of z/^'s if the thermalisation 
due to the mirror weak interactions is neglected. The unit along the vertical axis is MeV^. 

Figure 8: Evolution of the resonance momenta Pi/T and Pi/T for the example with Smi^^^^ = 
50 eV^ and Sm'^^^n — 1 qV'^- The bold lines on the right of the figure correspond to the 
^^fargc modes and the thin lines on the left of the figure correspond to the Sm^^^Yi modes. 
For the bold lines, the solid line, long-dashed line, short-dashed line, and dashed-dotted 
line correspond to Pi/T, P2/T, P^/T and P4/T respectively. For the thin lines, the solid 
line, long-dashed hne and dashed-dotted line correspond to pi/T,p2/T and ps/T ~ P4,/T 
respectively. 

Figure 9: Magnified version of the left-hand side of Figure 8. 

Figure 10: Evolution of the lepton numbers for the same example as Figs. 8 and 9. 

Figure 11: N^^cS versus ^w^smaii with Smf^j.^^ = 50 eV^. The case L^,^ < (which it turns out 
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implies L^^ > 0) has been considered. The dashed hne is the contribution SiN„^cS due to the 
effects of the Li,^ asymmetry while the dashed-dotted line is the contribution 62Nij^cS due to 
the change in the expansion rate. The solid line is the total contribution ^liVj, eff + 52-^i/,eff- 



Figure 12: 


Same 


as 


Figure 11 except 5m^^j.^^ — 200 eV^. 






Figure 13: 


Same 


as 


Figure 11 except Smi^^^^ — 800 eV^. 






Figure 14: 


Same 


as 


Figure 11 except L^^ > (and hence L^,^ 


<0) 


is considered. 


Figure 15: 


Same 


as 


Figure 12 except L,^^ > (and hence L^,^ 


<0) 


is considered. 


Figure 16: 


Same 


as 


Figure 13 except L^^ > (and hence L^^ 


<0) 


is considered. 
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